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ABSTRACT 

A  study  of  the  relationship  between  the  interferogram  F(x)  obtained 
in  Fourier  spectroscopy  and  the  corresponding  spectral  distribution  E (y) 
is  made,  it  is  found  that  the  spectral  distribution  is  the  Fourier  cosine 
transform  of  the  interferogram. 

Various  computational  techniques  for  performing  the  required  trans¬ 
formation  are  studied,  and  a  Real-Time  Fourier  Transform  Synthesizer 
is  described.  It  has  been  designed  to  be  compatible  with  a  lamellar 
grating  interferometer  but  could  be  used  just  as  well  with  a  well-compen¬ 
sated  Michels  on  interferometer.  The  system  used  to  obtain  the  point  of 
zero  path  difference  is  described  in  detail.  A  novel  approach  is  used  to 
obtain  the  cosine  functions  required  in  the  synthesis  of  the  cosine  transform. 
All  the  computations  except  the  summation  are  done  in  analog  form,  while 
the  summation  is  done  digitally.  The  result  of  the  computation,  the 
spectral  distribution,  is  available  as  soon  as  the  interferogram  is  obtained. 
Errors  in  the  technique  are  discussed  and  methods  of  compensating  for 
them  are  enumerated. 

Data  have  been  obtained  for  several  test  functions,  and  the  synthe  sired 
spectra  are  compared  with  the  theoretical  values.  Several  actual  intor- 
ferograms  have  been  obtained,  and  the  results  of  the  spectral  synthesis 
are  presented. 
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INTRODUCTION 

In  the  infrared,  there  are  two  reasons  for  using  an  interferometer. 
First,  an  interferometer  views  a  source  with  a  much  larger  entrance 
aperture  than  a  spectrometer  and  allows  more  energy  to  reach  the  detector. 
This  is  called  "aperture  gain. "  Second,  all  spectral  elements  in  a  band¬ 
width  determined  by  the  resolving  power  of  the  instrument  are  scanned 
simultaneously  during  the  full  time  of  the  observation.  If  T  is  the  obser¬ 
vation  time  per  spectral  element  and  N  is  the  number  of  elements,  it  is 
found  that  the  interferometer  makes  N  times  as  many  measurements  of 
each  spectral  element.  If  the  system  is  detector -noise  limited,  variations 
in  the  signal  due  to  random  detector  noise  vary  inversely  as  the  square 
root  of  the  observation  time  {  1]  and  this,  in  turn,  depends  on  the  number 
of  spectral  elements  being  observed.  This  fact,  the  gain  of  the  square 
root  of  N  in  the  signal-to-noise  ratio  of  the  recorded  signal,  is  known 
as  the  "Fellgett  advantage. "  Hence,  the  ratio  of  the  net  gain  of  an  inter¬ 
ferometer  over  a  spectrometer  is  the  aperture  gain  times  the  Fellgett 
advantage  [  2] . 

An  interferometer  may  be  used  as  a  spectrometer  by  measuring  the 
radiation  falling  on  the  central  fringe  of  a  two-beam  interferometer  as  a 
function  of  path  difference.  The  average  value  of  a  chart  recording  of  the 
radiation  as  a  function  of  path  difference  is  called  an  interferogram. 


The  spectral  distribution  may  be  obtained  by  calculating  the  Fourier 
transform  of  the  mterferogram.  A  digital  computer  is  generally  used  to 
perform  the  transformation*  The  use  of  a  digital  computer  has  a  dis¬ 
advantage  in  that  the  spectrum  is  not  available  for  immediate  analysis* 

This  research  is  undertaken  to  determine  the  feasibility  of  computing 
spectra  in  real  time  and  building  a  computer  to  perform  this  task. 
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TWO-BEAM  INTERFEROMETERS  AND  A 
DISCUSSION  OF  INTERFEROGRAMS 

Interferometers 

A  two-beam  interferometer  is  a  device  which  splits  an  incoming  wave- 
front  into  two  parts  by  amplitude  or  wavefront  division.  It  introduces  a 
delay  in  the  path  of  one  part  and  then  recombines  the  two  part  s.  The 
Michelson  interferometer  is  an  example  of  the  type  which  operates  by 
amplitude  division,  see  Figure  1  .  The  beamsplitter,  providing  it  is 
non -absorbing  and  has  a  reflectance  of  0.  5,  divides  the  light  into  two 
equal  parts  and  recombines  them  after  a  relative  delay  of  or  thus  causing 
them  to  interfere.  The  delay  is  produced  by  introducing  different  optical 
paths  for  the  different  parts  of  the  divided  beam.  One  mirror  is  fixed 
while  the  other  can  be  displaced  parallel  to  itself  with  a  constant  velocity 

i 

v.  If  M  is  shifted  through  distance  d  to  position  M_  ,  the  time  delay 
£  » 

associated  with  this  path  is  increased  by  2  d/c  where  c  is  the  speed  of 
light.  The  geometrical  path  is  increased  by  2d.  The  optical  path  is 
increased  by  2md  where  m  is  the  index  of  refraction.  A  compensator  is 
included  in  one  path  to  equalize  the  paths  in  material  other  than  air.  The 
Michelson  interferometer  exhibits  circular  symmetry,  and  the  detector 
is  placed  in  the  central  fringe  of  the  interference  pattern. 
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An  interferometer  which  makes  use  of  wavefront  division  is  the 
lamellar  grating.  Figure  33  show*  a  system  in  which  the  lamellar 
grating  is  used,  and  Figure  17  shows  a  cross  section  through  the  facets. 
The  odd-numbered  facets  are  fixed,  and  the  even-numbered  ones  are 
mounted  in  such  a  manner  that  they  are  movable  over  a  total  distance  of 
one  inch.  In  this  particular  case,  both  sets  of  facets  are  made  of  steel 
although  other  materials,  such  as  glass,  have  been  used.  The  grating  is 
6  inches  square,  and  each  facet  is  6  inches  long  and  1/16  inch  wide.  The 
faces  of  the  facets  have  a  deposited  gold  coating  to  increase  the  reflectance 
in  the  infrared.  Each  facet  acts  as  a  narrow  mirror,  and  the  path  differ¬ 
ence  is  introduced  by  changing  the  groove  depth.  Interference  takes  place 
in  a  plane  located  at  the  faces  of  the  foremost  facets.  The  detector  is 
placed  in  the  central  grating  order. 

Lamellar  grating  interferometers  are  used  in  the  infrared  region  of 
the  spectrum.  The  one  used  in  this  research  is  designed  for  wavelengths 
longer  than  10  microns.  They  are  not  used  in  the  visible  region  because 
it  is  impossible  to  machine  and  polish  the  surface  of  the  facets  to  the 
required  tolerances. 

The  illumination  received  by  the  detector  is  modulated  by  translating 
the  movable  section  of  the  grating  or  one  mirror  of  the  Michelson  interfero¬ 
meter.  When  this  is  done  linearly  with  time,  the  detector  response  is 


modulated  at  a  frequency 
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=  2w  v 
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where  v  is  the  velocity  of  the  mirror,  or  moving  facets,  and  vq  is  the 

wavenumber.  Also,  v  ■  1/h  ,  where  ia  the  wavelength. 

o  o  o 

For  monochromatic  illumination,  the  detector  response  ia  sinusoidal 
and  would  continue  indefinitely.  When  several  monochromatic  components 
are  present,  the  detector  response  is  the  sum  of  a  series  of  sinusoidal 
functions,  one  for  each  of  the  components.  Thus,  the  output  of  an  inter¬ 
ferometer,  called  an  "interferogram,  "  is  synthesized  by  means  of  the 
summation  of  trigonometric  functions  with  amplitudes  determined  by  the 
spectrum  of  the  incident  illumination.  By  computing  the  Fourier  transf: -  vm 
of  the  interferogram,  it  ia  possible  to  make  a  harmonic  analysis  of  it  and 
thence  to  determine  the  intensities  of  the  frequency  components.  Thus, 
the  spectrum  may  be  obtained. 

Derivation  of  equations 

In  the  infrared,  one  is  dealing  with  frequencies  on  the  order  of  10*^- 
15 

10  cycles  per  second.  Present  detectors  cannot  respond  to  frequencies 

7 

higher  than  about  10  cps.  When  higher  frequencies  fall  upon  these 
detectors,  they  respond  to  the  time  average  intensity  of  the  illumination. 
Time  dependence  will  be  ignored  in  the  following  work,  although  it  can  be 
shown  that  the  time  origin  is  arbitrary. 
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The  following  derivation  it  applicable  to  the  case  of  the  Michelson 
interferometer  or  the  lamellar  grating.  Consider  an  electromagnetic 
wave  of  single  frequency  v  , 

-i2irv  6 

v(6}  » voe  °  (l ) 


where 

v  is  in  wavenumbers  and  is  equal  to  1  f\  , 
o  o 

is  the  wavelength  in  centimeters, 
vq  is  the  amplitude  of  the  wave,  and 
6  is  the  distance  from  some  arbitrary  origin. 

The  beamsplitter,  or  the  grating,  divides  the  wave  into  two  parts 
Vj{6)  and  v2<6)  such  that  av^  =  For  the  Michelson,  "a"  may  be  a 
constant,  a  comp! ex  quantity,  or  a  frequency  dependent,  as  determined 
by  the  characteristics  of  the  beamsplitter.  It  is  nearly  unity  for  the 
lamellar  grating. 

Let  the  difference  in  optical  path  length  between  v  (6)  and  v  (6),  na 
measured  from  the  beamsplitter  or  the  front  face  of  the  grating,  be  x. 
Then 


-iZtrv  5 

Vj(6)  =  vte  ° 


s 

w 


! 


4 

i 

I 

V 


(2) 
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and 

-i2irv  (5  +  it) 

v2(4)=avj*  °  (3) 

At  the  recombination  point,  the  cum  of  these  is 

-i2irv  6  -i2wv  x 

v(6)  »  +  v2(6)  *  VjC  °  (i  +  ae  °  )  (4) 

This  equation  demonstrates  that  the  origin  of  6  is  arbitrary  and  may  be 
chosen  equal  to  aero. 

At  any  point  in  space,  the  intensity  is  given  by 

I  »  v(6)  v*(6)  (5) 


where  v*{6)  represents  the  complex  conjugate  of  v(6).  Therefore, 

i2trv  x  -i2irv  x 

Is  V1  [1  +a  +a(e  +  e  °  )]  (6) 

For  an  ideal  beamsplitter,  v  =  v  ,  so  "a"  s  1  in  this  case  also.  Since 

*  c* 

it  is  only  a  function  of  x,  equation  {6}  may  be  rewritten  as 
I(x)  s2vj[I  +  cos  2wvox] 


(7) 


o 


and 
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v(v)  s  £(v)  s  intensity  spectral  distribution. 


With  these  substitutions,  equation  (11)  may  be  rewritten  as 


E(x)  *  f  E(v)[l  +  cob  2trvx]  dv> 


(12) 


I2(x)  =  f  E(v)dv  +  ^  £{v)  cos  2ir  xdx. 


(13) 


For  x  =  0 


os 

E(0)  =  2  f  E(v)dv  (14) 

and  lor  large  x  *  d 

00 

E(d)  a  |E(0)  +  f  E(v)  cos  2TTVd  dv  (15) 

“00 

The  integral  will  vanish  unless  E(v)  represents  strictly  monochromatic 
radiation.  In  practice,  the  value  of  d  required  for  the  integral  to  vanish 
is  determined  by  the  wavelength  range  of  the  incident  illumination.  The. 
broader  the  bandwidth,  the  smaller  the  value  of  d  required  to  yield  E(u)  - 
iE(0). 

The  difference  between  E(x)  and  |e(0)  is  analysed  and  is  called  the 
"interferogram  function"  F(x). 


QO 


F{x)  *  E(x)  -  |E(0)  mj  E(v)  cos  Zirvxdv 


This  function  is  obtained  from  a  strip  chart  recording  of  the  total 
illumination  falling  upon  the  detector  of  the  two -beam  interferometer 
as  the  path  difference  is  varied.  The  value  of  |E(0)  is  determined  from 
the  asymptotic  value  of  such  a  record.  The  total  illumination  approaches 
this  value  as  the  path  difference  increases. 

From  Fourier  theory 


E(v)  s  t  F(x)  cos  2wxdx 


Therefore,  if  we  can  specify  F(x),  we  can  uniquely  determine  E(v). 


Methods  of  coi 


spectra  from  interferograms 


For  simple  interferograms  which  can  be  expressed  analytically,  the 
computation  can  proceed  directly.  Two  such  cases  are  discussed. 

1)  Assume  the  interferogram  can  be  expressed  as 

_  2.  2 

F(x)  *  e  X  a  cos  bx  (18> 


The  spectrum  is 


00  2  2 

E(v)  s  jTf(x)  cosZirvxdx  *  ^  COB  k*  cos2vvxdx  (19) 


00 


1? 


This  may  be  solved  with  the  aid  of  integral  tables  and,  by  neglecting 
constants,  the  result  is 


E(v)  s  e 


2.  ,  ,2 

•a  tv-b)  /4 


<20) 


This  represents  a  line  with  a  Gaussian  frequency  distribution  centered 
at  wavenumber  b. 

2)  If  the  interferogram  can  be  expressed  as 


F(x)  s  e 


2.  2 

-x  /  a 


cos  fx  cos  bx 


where 


b  »  f ,  af  »  1 


(21) 


then 


00  2  ,2 

-x  /a 


E(v)  a  J'  F(x)cos2trvx  dx  m  J  e”X  a  roc  C:  cos  bx  ecc  2irvxdx  (22) 


This  may  also  be  solved  with  the  aid  of  integral  tables  and,  by  neglecting 
constants,  the  result  is  [3] 


E(p) 


-a2(v-b-f)2/4  ,  -a2(v-b+f)2/4 

e  '  ’  +  e  ' 


(23) 
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This  represents  a  pair  of  lines  with  Gaussian  distributions  centered  at 
wavenumbers  b-£  and  b+£.  The  interferograms  and  the  corresponding 
spectral  distributions  are  shown  in  Figure  2. 

The  interferogram  function 

00 

F(x)  a  E(X)  -  £E(0}  =  f  E{v)  cos2irvxdv  (24) 

-00 

may  be  rewritten  in  terms  of  electrical  engineering  terminology  in  terms 
of  w,  where  u  a  2'rf  and  is  the  angular  frequency  of  the  detector  output. 

F urther,  a  2vv  where  v  a  the  rate  of  change  of  groove  depth  or  mirror 
position.  Let  cp(w)  a  ~  E(v)  and  G(t)  a  F(x)  [4],  Then 

00 

G(t)  ^  q>(w)  cosutdu  (25) 

”te 

where  cp(u)  is  the  frequency  spectrum  of  the  detector  response  and  E^(v) 
is  the  spectrum  of  the  radiation  incident  on  the  detector.  The  frequency 
spectrum  is 


<p(u}  a  f  G(t)  coswtdt 


(26) 


and  the  problem  is  reduced  to  that  of  finding  qp(u)  if  given  a  pi r'c  of  G(t). 
Since  G(t)  cannot  generally  be  expressed  analytically,  some  sort  of 
analog  or  digital  computer  is  used  to  obtain  the  Fourier  transform. 


asMBKiea  -  v 


m  stum* 


•  15 

Analog  computers  have  been  used  to  compute  Fourier  coefficient*, 
ever  since  the  Hendrici  Analyzer  was  constructed  in  the  early  part  of 
this  century  £  5]  *  Many  other  computers  were  constructed,  but  most  cf 
them  were  mechanical,  slow,  and  of  limited  accuracy  [6,  7,  8,  9] .  Sever?,! 
optical  devices  have  been  investigated,  but  they  do  not  lend  themselves  to 
the  real-time  calculation  of  spectra  from  interferograms  [  10,  11].  One 
of  these  optical  schemes  is  discussed  in  the  Appendix.  Firle  and  Staake 
recorded  data  on  film  or  tape  and  played  it  back  through  a  wave  analyzer 
£12, 13]. 

Neither  the  conventional  sweep  analyser  nor  the  bank-of -filters 
analyzer  will  work  in  the  very  low  frequency  range  of  the  interferogram. 
When  these  units  are  used,  the  general  technique  is  to  record  the  interfer¬ 
ogram  on  tape  and  play  it  back  at  a  higher  speed,  thus  increasing  all  the 
frequencies  by  the  speed-up  ratio.  Analyses  by  this  method  are  either 
quite  slow  or  have  limited  resolution. 

The  IFTAC  ( 14]  records  the  interferogram  on  a  tape  loop  and 
actually  computes  one  component  each  time  the  tape  loop  is  played  back. 

It  is  designed  to  compute  the  spectrum  and  to  give  1000  output  points. 

The  SIMORAMIC  Analyzer  operates  much  faster  than  a  conventional 
sweep  analyzer  of  comparable  bandwidth  and  resolution  [15].  It  hac 
been  successfully  used  in  the  range  of  0-8  kc  to  study  atmospheric 
whistlers  [  16].  Until  the  advent  of  high  speed  digital  computers, 
interferometry  fell  into  disuse,  partially  because  of  the  laborious  task 


of  analysing  the  data  and  the  inadequate  computers  and  analyzers  then 
available. 

Only  discrete  values  of  G(t)  can  be  utilized  by  a  digital  computer  or 
the  r eal -time  Fourier  transform  synthesizer.  Therefore,  to  compute 
cp(w),  it  is  necessary  to  sample  G{t).  The  samples  must  be  spaced  close 
enough  to  obtain  all  the  information,  but  not  so  close  as  to  add  complicate” 
without  a  gain  in  information.  The  choice  of  the  sampling  points  is 
governed  by  Shannon's  sampling  theorem  [  17]  •  The  sampling  theorem 
in  the  time  domain  is: 

If  a  band  limited  function  G(t)  contains  no  frequencies  higher 
than  u0  radians /second,  it  is  completely  specified  by  giving 
its  ordinates  at  a  series  of  points  spaced  ir/w0  seconds  apart, 
the  series  extending  throughout  the  time  domain. 

Thus,  if  the  maximum  frequency  of  the  detector  response  is 

uq/4v  v/X0»  it  is  sufficient  to  sample  the  interferogram  at  intervals 

Ah  a  h  /4  =  l/4v  seconds.  The  range  of  values  which  v  can  assume  m?y 
o  o 

be  selected  by  optical  filters. 

The  sampling  theorem  in  the  frequency  domain  is  important  since  it 
tells  us  at  what  interval  in  w  rad/sec.  (or  vcm"*)  we  must  compute  the 
spectral  function  <?(«}  (or  Ej(v))  to  completely  specify  it.  The  theorem 
is  stated  as  follows: 

If  q>(w)  represents  the  spectrum  of  a  function  G(t)  which  is  zero 
everywhere  except  in  the  interval  (tj  <  T,  then  q>( u>)  is  exactly 
determined  for  all  values  of  «  by  giving  its  values  at  series  of 
points  1/T  cycles /second  apart  in  frequency,  the  series 
extending  throughout  the  frequency  domain  [18], 
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No  information  is  lost  if  is  specified  at  points  separated  by  tr/T 
radians  or  1/4D  in  wavenumbers  where  D  is  the  maximum  groove  depth 
or  mirror  displacement. 

The  sampling  theorems  tell  us  the  proper  sampling  interval  at  which 
we  do  not  lose  any  information.  It  is  possible  to  sample  at  faster  rates, 
and  it  is  generally  done  in  practice,  as  signal -to -noise  ratios  are  improve  , 
If  h  is  the  sampling  interval,  the  sampled  spectra  repeat  at  intervals  of 
i/h*  Since  noise  is  present  in  all  practical  cases,  the  sampling  rate  is 
increased  on  the  order  of  5  to  prevent  overlapping  of  noise  in  the  repeating 
spectra  [19]. 

When  digital  computers  are  used  to  calculate  the  Fourier  transform, 
they  are  generally  programmed  to  analyze  the  interferogram  F(x)  to 
obtain  the  spectrum  E(v).  That  is  to  say,  every  sample  point  is  used  to 
obtain  information  about  each  component  of  the  spectrum.  In  the 
synthesis  technique,  used  in  the  real-time  Fourier  synthesizer,  each 
sample  point  gives  some  information  about  all  the  speciral  components. 

Assume  the  interferogram  G{t)  is  sampled  at  the  proper  intervals 
from  -T  <  t  <  T.  This  represents  a  series  of  impulses  with  amplitudes 
determined  by  the  value  of  the  interferogram  at  the  sampling  point,  and 
the  spacing  is  dependent  on  the  sampling  interval,.  Figure  3  .  The  in?-c?- 
ferogram  G[t)  is  then  represented  approximately  by 
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G(t)  a  g{ 0)6{0)  +  gfh)d(t-h)  +  g(-h)5(t+h)  +  g(2h)6(t-2h)  +  g(~2h)6(t+2h) 

+  ....+  g(nh)6(t-hn)  +  g(-nh)<5(t+nh)  (27) 

The  Fourier  transform  of  a  pair  of  6  functions  symmetrically  placed 
about  the  origin  is  a  cosine  function  whose  frequency  depends  on  the 
distance  of  the  pair  of  6  functions  from  the  origin*  Therefore,  by  super¬ 
position,  the  Fourier  transform  of  G(t)  is 

<p(«)  *  g(0)  +  g(h)  coso^t  +  g(2h)  cos«2ht  +.  ,  +g(nh)  coso^t 

n=N 

<p(w)  *  g(0)  +  r  g(nh)  cosw^t  (28) 

n=0 

This  equation  states  that  the  spectrum  may  be  synthesized  by  the 
summation  of  a  series  of  cosine  functions  with  amplitudes  dependent  upon 
the  amplitudes  of  the  sample  and  frequencies  dependent  upon  the  distance 
of  the  sample  point  from  t  =  0. 

Scanning  function 

In  practice  it  is  not  possible  to  sample  the  interferogram  throughout 
the  entire  time  domain,  and  a  truncation  of  the  interferogram  results 
The  Fourier  series  representation  of  the  spectrum  estimate  in  exponentlol 


form  is 
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n«N  .  , 

<PN(u)  h  h  S  (29) 

n«-N 


This  converses  to  the  actual  spectrum  as  N“* oo,  but  since  the 
interferogram  is  truncated,  only  a  finite  number  of  terms  are  available. 
This  will  cause  the  upectrum  estimate  to  deviate  from  the  true  spectrum. 
The  computed  spectrum  will  be  the  convolution  of  the  true  spectrum  with  a 
function  called  the  scanning  function  which  will  now  be  derived. 

In  equation  {29}  let 


"o 

G{nb)  »  3“  f  <?{«)  ejnhi  d* 

**C£ 

O 

Then,  equation  {29}  becomes 


°  N 

<Pa<w)“h2~  J  <PC«)  Sejn(<_w)hd*  (31' 

-u  ^ 

o 

Since 


is  e  «  i+cos{« -o>)h  +  cos2{«  «%j)h+. . +cosN{« -w)n  {32} 

n*-N 


the  right-hand  part  of  equation  (32)  is  equal  to  [20] 


sin(2N+l)  («-w)h 
_  2 

sin  (^r-)h 


(i,  («  -W)  «8 
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Substituting  equation  (33)  into  equation  (31)  yields 


<Pn(w)  = 


Jl 

2n 


-  u>)dt 


u»o 

-  f  <p(«]p(«-co)tU 


(34) 


For  large  N,  the  factor  1  in  equation  (33)  is  insignificant  and  can  be 
dropped.  Also,  since  the  interferogram  is  even,  only  N  samples  are 
independent,  so  equation  (33)  can  be  rewritten 


<«> 

This  is  called  the  scanning  function.  In  electrical  engineering  terminology, 
it  would  be  called  the  impulse  response.  From  equation  (34),  it  is 
observed  that  the  spectrum  estimate  is  the  convolution  of  the  true 
spectrum  with  the  scanning  function  as  the  integration  is  carried  from 
0  to  (jjq,  the  fundamental  range.  The  estimate  q>n(«)  approaches  cp(a>) 
more  closely  as  N  increases.  The  first  aero  of  ;j .(«  -w)  occurs  when 


A  co  * 


2  IT 

Nh 


Thus,  for  a  monochromatic  source,  it  is  found  that  the  computed  sper^ra'. 
line  has  a  width  that  is  inversely  proportional  to  N. 

For  a  pair  of  closely  spaced  lines,  the  Rayleigh  criterion  gives  a 


resolving  power  [21] 
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It  is  now  useful  to  see  how  the  spectrum  estimate  approximates  the 
true  spectrum  in  the  region  of  a  discontinuity.  In  Figure  4  ,  the  true 
spectrum  is  assumed  to  be  a  square  pulse  and  the  scanning  function  is  as 
shown.  The  computed  spectrum  estimates  are  given  for  two  values  of  N. 

The  scanning  function  has  a  large  peak  and  several  smaller  peaks  on 
either  side  which  decrease  in  amplitude.  The  negative  peaks  adjacent  to 
the  main  peak  have  amplitudes  slightly  greater  than  one -fifth  the  value  cf 
the  main  peak,  and  they  do  not  become  smaller  no  matter  how  large  N 
becomes. 

The  sine  integral  of  x  is  shown  in  Figure  5  .  It  is  seen  that  the 

integral  never  assumes  its  final  value  but  keeps  continuing  to  approach  it. 

The  peaks  occur  at  x  =  it,  5v,  5ir.  . .  etc. ,  and  the  valleys  occur  at  2tr,  4ir, 

etc.  The  separation  «.!  the  peaks  in  x  is  2r.  This  gives  a  separation  in  w 
2 

of  Aw  3  * 

In  the  case  where  the  smaller  peaks  become  compressed  into 

a  single  vertical  line  at  a  discontinuity,  but  even  then  the  Fourier  ser?r^ 
still  yields  a  9%  overawing  compared  to  the  value  of  the  discontinuity. 
Thus,  to  reduce  this  error  due  to  the  Gibbs  phenomena  at  a  discontinuity, 
such  as  a  distinct  spectral  line,  a  technique  called  "apodizatior’1  i*»  u"cd. 
This  can  be  seen  more  clearly  by  the  following  derivation. 
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Figure  4 


The  true  spectrum,  in  terms  of  wavenumbers,  i»  given  by 

i 

00 

E(v)  a  f  F(x)  cos2irvx  dx  {3”) 

j 

kJ 

”09 

The  maximum  path  difference  is  finite,  and  this  causes  truncation  of  the 
interferogram  function  F(x),  Equation  (37)  may  now  be  expressed  as 

CO 

En(v)  =  j'  F(x)  C(x)  cos2itvx  dx  (38) 

•oo 

where 


C(xj  a  1  for  -EKx<+D 
C(x)  =  0  for  x<  -D  and  x>  +D 


and 


D  s  maximum  path  difference. 

The  cosine  transform  of  the  product  of  two  even  functions  is  the  cotivoV-J  ,  .■ 
of  the  cosine  transforms  of  each  of  the  two  functions.  Therefore, 

E^(v)  s  [cosine  transform  of  F(x)J  *  [cosine  transform  of  C(. ..1 . 

The  cosine  transform  of  F(x)  is  the  true  spectrum,  while  the  cocin<> 
transform  of  C(x)  is  the  scanning  function.  The  scanning  function  hr.'-  e 


'I  I*- 
*  w 


c’.lmcnsion  of  1/length  when  E(v)  is  computed  and  l/tlme  .-'hen  -p'w)  ic 
computed.  Now  define  the  cosine  transform  of  C(x)  =  n{«  -  ■  ),  This  gives 
the  same  result  as  previously  obtained. 

If  the  interferogram  function  is  multiplied  by  a  function  A{x)  other  than 
C(x),  equation  (33)  becomes 


TO 


F{x)  A(x)  costrvxdx 


(40) 


A(x)  is  called  an  apodizing  function  and  is  chosen  to  possess  the  property 
of  having  low-side  lobes*  A  commonly  used  function  is  the  triangle  sh 
in  Figure  6  . 


Aj(x)  *  1  -  x/d  for  -D<  x  <  D 

A(j-)  »  0  for  x  <  -D  and  x  >  D 


The  Fourier  transform  of  this  function  is  of  the  form 


p.{<  -x)  *  D  (sin 


rwD  2 
i»D' 


It  has  no  negative  ripples  and  much  smaller  side  lobes  than  the  sine 
function  obtained  for  square  pulse  apodieation. 

Another  useful  apodizing  function  is  of  the  form 


A2(x)  a  [1«x2/D2]2 
a  0 


for  «D  <  x  <  D 
for  x  <  -D  and  x  >  D 


APODIZING  FUNCTIONS  AND  THEIR  CORRESPONDING  SCANNING  FUNCTIONS 
4c<x) 


For  this  case,  the  Fourier  transform  is  of  the  form  [23] 


„c/2 

!±(«-x}*{2  *Dv)  Js/zttrrDv) 

Figure  6  shows  this  function  and  the  effective  scanning  function. 
Resolving  power 

The  terms  "spectral  element"  and  "resolution"  have  been  used.  Let 
us  now  examine  their  meaning.  The  resolution  or  resolving  power  of  a 
system  is  a  measure  of  how  close  two  adjacent  lines  or  sources  can  be 
before  they  are  no  longer  sufficiently  separated  so  as  to  be  identified. 
Assume  a  pair  of  lines  of  arbitrary  intensity  and  separated  by  Av  cm 
The  lines  will  be  considered  to  be  resolved  if  the  difference  between  tv,e 
intensity  of  the  weaker  line  and  the  intensity  where  it  overlaps  the  stro-v 
line  is  at  least  twenty  per  cent  of  the  intensity  of  the  weaker  line.  For 
two  lines  which  can  be  considered  as  delta  functions,  the  overlap  is 
produced  by  the  scanning  function  of  a  spectrometer,  or  the  scanning 
function  due  to  apodization  for  the  case  of  an  interferometer.  For 
triangular  apodization,  this  is  equivalent  to  the  Rayleigh  criterion.  This 
means  that  the  scanning  function  for  a  line  at  v  cm  *  has  a  maximum 
while  the  scanning  functions  for  lines  at  v  +  dy  cm  *  have  zeros  at  v  cm 
Figure  7  . 


For  a  triangular  apodizing  function  of  length  D,  its  Fourier  trail.; '’--im 

2 

(the  scanning  function)  is  sine  (uvD),  Two  adjacent  lines  scanned  by  tb.‘> 
function  overlap  at  a  point  midway  between  them,  where  the  value  is 
approximately  0.  8, 

Figure  6  shows  several  apodizing  fun  dons  and  their  corresponding 
scanning  functions.  The  first  zero  of  the  scanning  function  for  rectangular 
truncation  (no  apodization)  occurs  at  half  the  frequency  of  the  first  zero  of 
the  scanning  function  for  triangular  apodization.  That  is,  the  resolving 
power  is  increased  by  a  factor  of  two  when  no  apodization  is  used.  The 
spectrum  in  this  case  will  have  spurious  oscillations  due  to  the  oscillat’nr  <3 
of  the  scanning  function.  The  Fourier  transform  of  a  function  gets 
narrower  as  the  function  gets  wider,  and  one  can  see  immediately  that  tb° 
limit  of  resolution  varies  inversely  with  the  maximum  path  difference. 

The  exact  variation  is  determined  by  the  particular  case.  The  limit  of 
resolution  will  be  callod  a  spectral  element.  That  is,  it  is  the  emalleet 
value  of  Am  by  which  two  lines  can  differ  and  still  be  resolved. 


REAL-TIME  FOURIER  COSINE  TRANSFORM  SliJTIJEr  — R 
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Principles  of  Operation 


The  interferogram  is  represented  by 


v 

F(x>*  f 


M 

E(v)  cos2wx  dv 


(<•1) 


v 

m 

where  v  and  ,  are  the  minimum  and  maximum  wavenumbers  in  the 
m  M 

spectral  distribution*  It  is  evident  from  equation  (41)  and  from  the 
discussion  in  the  preceding  chapter,  that  the  interferometer  has  syniher  Vr 
the  inteiferogram  by  a  superposition  of  trigonometric  functions  with 
amplitudes  determined  by  the  spectral  distribution.  Each  point  on  the 
interferogram  may  be  represented  by  a  finite  polynomial,  the  number  at 
terms  being  limited  by  the  upper  frequency  limit  of  the  radiation*  This 
polynomial  is  a  partial  sum  of  a  Fourier  series  [24]. 

The  spectral  distribution  is  related  to  the  interferogram  by 


00 

E(v)  =  ^  F(x)  cos2irvxdx  (42) 

"V> 

It  is  reasonable  then  to  obtain  the  spectral  distribution  by  the  inverse 
process;  i.e*,  by  a  superposition  of  trigonometric  functions  with 
amplitudes  determined  by  the  interferogram*  A  synthesis  of  the  spectre! 
function  is  thus  obtained.  This  synthesis  is  the  Fourier  series  approxi¬ 
mation  to  the  spectrum. 


The  interferogram  function  F(x)  is  sampled  v t  ?  .;er.‘*is  jerminec 
by  the  sampling  theorem*  F{x)  may  then  be  expressed  as  F*(x),  where 

ee 

F »(x)  »  E  A(x-nh)  <43} 

ns-oo 

« 

The  expression  E  A(x-nh)  is  a  Dirac  comb  of  periodicity  h.  It  has  the 
property  that  it  picks  up  values  of  F(x)  at  equal  intervals  h  of  x,  where  h 
is  die  spacing  between  sample  points.  Its  Fourier  transform  is  also  p. 

Dirac  comb,  but  of  periodicity  1  /'h*  Further,  F(x^)  is  the  value  of  F(x- 
at  sample  point  nh,  see  Figure  8  * 

Equation  (42}  may  now  be  written  as  a  summation  instead  of  an  integral. 

CO 

E(v)  *  E  F(xwK)  cos2irvxwh  f 

n«-v»  1 

Since  the  interferogram  has  been  sampled,  the  synthesized  spectrum  i*? 
the  convolution  of  the  true  spectrum  with  the  Dirac  comb;  thus,  the 
synthesized  spectrum  repeats  at  intervals  of  1  /h. 

Equation  (44)  requires  an  infinite  number  of  samples,  which  is 
clearly  a  practical  impossibility.  Path  differences  are  governed  by  the 
physical  limitations  of  the  instrument  and  only  a  finite  number  of  aar.  Aoo 
are  obtained.  This  limits  the  resolution  of  the  instrument,  but  it  mav  st-’l 
be  very  large  because  the  path  difference  can  be  made  much  larger  th»r> 


the  ong  wavelength*  The  in.  ^rfer  -g-am  ij  \n  »v  .  iv^  "ti  .  'o  oj  v 
t  .ote  sample  points  from  n  ■  0  to  N  need  be  used,  »vhero  F  ip  '‘he  "net 
sample  point  and  is  governed  by  the  maximum  path  difference.  With  '*;ese 
stipulations*  equation  (44)  may  be  written 

N 

E(v)  «  2  **{*31*)  co^irvx^  (45) 

nsO  “  4 

The  relationship  between  the  functions  E(v)  and  qp{w)  and  between  F(x)  and 
G(t)  have  been  discussed  previously.  Equation  (45)  may  be  rewritten  ir 
terms  of  these  parameters  as 

N 

<p(t*>)  *  E  G(t.  )  cosu^t  (46) 

n»0  n 

where  cp(a)  is  the  spectrum  of  the  detector  response*  and  G(t  .  }  is  the 

nil 

amplitude  of  the  detector  response  at  the  sample  point  nh.  The  angular- 
frequency  of  the  cosine  functions  is  given  by 


u  =  2i rf 
n  n 


where  is  the  bandwidth  of  the  incident  radiation  and  the  running  va- 
t  is  equivalent  to  the  parameter  v.  Each  sample  point  of  the  interferogi  •- 
is  used  to  generate  a  cosine  function  with  an  amplitude  determined  by  th 
sample  value  and  a  frequency  determined  by  the  path  difference. 


3' 


Equation  (46)  forma  the  basis  of  the  real-time  i  ouricr  Ermthesizor. 

The  sample  value  G(t^)  is  held  while  it  is  multiplied  by  a  cosine  function 
of  the  proper  frequency.  Each  succeeding  product  is  put  into  the  memory 
and  is  added  to  the  preceding  data.  A  summation  of  G(t^K)  cosu^t  terms 
is  thus  obtained. 

A  photographic  summation  technique  was  investigated.  It  has  some 
severe  disadvantages  which  limit  its  usefulness.  A  discussion  of  the 
technique  and  the  results  of  the  experiments  are  given  in  the  Appendix. 

As  a  result  of  this  investigation,  the  technique  was  abandoned  in  favor 
of  a  f tally  electronic  system. 

A  photograph  of  the  real-time  Fourier  transform  synthesizer  is 
shown  in  Figure  9  and  an  elementary  block  diagram  is  presented  in 
Figure  10  .  Four  signals  are  required  to  synthesize  the  spectral  distri¬ 
bution.  Three  of  these  are  inputs  to  the  transform  synthesizer.  They 
are:  the  output  of  the  lamellar  grating,  the  reference  signal  from  the 
Michels  on,  and  the  initiating  pulse  from  the  zero  path  difference  det-rctor. 
The  cosine  functions  are  generated  within  the  instrument. 

The  relationship  between  the  lamellar  grating  and  the  Michelson  ?3 
shown  schematically  in  Figure  33  ,  while  Figure  11  is  a  photograph  o' 
the  two  interferometers  and  the  driving  system.  The  movable  facets  of  t.  ? 
lamellar  grating  are  designated  by  the  letter  "D"  while  the  fixed  facets 
are  designated  "E.  "  Mirror  of  the  Michelson,  is  mounted  on  the 
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REAL  TIME  FOURIER  SYNTHESIZER 


FRONT 


REAR 


Figure  9 


REAL  TIME  FOURIER  SYNTHESIZER 
BASIC  BLOCK  DIAGRAM 


Figure  10 
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LAMELLAR  GRATING  AND  MICHELSON  INTERFEROMETERS 


D  =  MOVABLE  FACETS 
E  =  FIXED  FACETS 
Mi  -  FIXED  MIRROR 
M2  =  MOVABLE  MIRROR 
B  -  BEAMSPLITTER 
C  =  CONPENSATOR 


HG  =  MERCURY  LAMP 
F  =  5461 X  FILTER 
Mo  =  MOTOR 

MS  =  MICROMETER  SCREW 
Ph  =  PHOTOMULTIPLIER 


Figure  11 
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baseplate  while  is  mounted  on  the  slide  which  carries  the  movable 
facets.  Items  "B*1  aad  "C"  are  the  beamsplitter  and  compensator, 
respectively.  A  mercury  lamp  and  filter  are  the  source  of  a  5461  k 
radiation  which  is  detected  by  a  photomultiplier.  A  1  / 1 0  rpm  synchronous 
motor  and  micrometer  screw  are  used  to  drive  the  interferometers. 

The  output  of  the  photomultiplier  varies  sinusoidally  as  a  function  of 
mirror  displacement.  Each  cycle  indicates  a  change  in  displacement  of 
\/2„  or  0. 27305  microns.  This  signal  is  used  to  measure  displacements, 
and  its  accuracy  is  independent  of  any  imperfections  in  the  driving  system. 
This  reference  signal  is  used  to  control  the  sampling  intervals  and  the 
frequency  of  the  cosine  functions. 

In  order  to  synthesise  the  Fourier  cosine  transform  of  the  inter* 
ferogram,  it  is  necessary  to  know  when  the  path  difference  between  the 
two  beams  is  aero.  This  information  is  supplied  by  the  aero  path  difference 
detector.  A  novel  approach  is  used  to  generate  a  pulse  when  the  groove 
depth  of  the  lamellar  grating  passes  through  zero.  The  approach  and  the 
results  are  explained  more  fully  in  a  section  devoted  to  this  particular 
circuit.  The  pulse  initiates  the  automatic  operation  of  the  transform 
synthesiser. 

There  are  four  main  signal  paths  in  the  instrument.  Two  of  these 
initially  end  at  an  open  relay  while  the  others  are  complete  but  inactive. 

The  aero  path  difference  detector  generates  a  pulse  which  closes  the  relit*/. 
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This  complete  a  the  paths  taken  by  the  outputs  of  the  two  interferometers 
and  initiate*  the  sampling  of  the  iaterferogram. 

With  reference  to  Figure  10  ,  the  first  path  begins  at  the  output  of 
the  lamellar  grating,  extends  through  the  sample  and  hold  circuit  and 
ends  in  the  multiplier.  The  sampler  takes  a  sample  of  the  interferometer 
signal  and  holds  it  during  the  interval  between  successive  samples. 

During  this  period,  the  sample  value  is  multiplied  with  the  output  of  the 
second  path  in  an  analog  multiplier. 

The  proper  cosine  functions  are  generated  in  the  second  signal  path. 
This  begins  at  the  frequency  synthesizer,  extends  through  the  gate  circuit 
and  ends  in  the  multiplier.  The  frequency  synthesizer  generates  a  sinu¬ 
soidal  signal  of  known  frequency.  It  is  remotely  programmed  and 
automatically  changes  frequency  according  to  a  pre-set  program.  The 
gate  opens  at  the  proper  instant  and  cosine  functions  enter  the  multiplier. 

The  first  two  signal  paths  are  joined  in  the  analog  multiplier  to  form 
the  third.  The  sample  values  G(t^)  are  multiplied  by  the  proper  cosine 

functions  to  produce  a  series  of  terms  of  the  form  G(t  ,  )  cosw  t.  These 

nn  n 

terms  follow  the  third  path  to  the  summation  unit. 

Figure  12  is  a  diagram  which  shows  the  waveforms  in  these  fiist 
three  signal  paths.  The  input  and  output  of  the  sample  and  hold  circuit 
are  shown.  The  two  inputs  to  the  multiplier  and  the  output  of  it  are  also 
presented.  Figure  13  shows  the  result  of  the  summation  of  these  first 
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few  terms.  The  synthesized  spectrum,  represented  by  the  summation, 
is  seen  to  repent,  as  would  be  expected  from  a  consideration  of  the 
sampling  process. 

The  fourth  signal  path  in  i  igure  10  is  the  one  taken  by  the  reference 
signal.  It  extends  from  the  detector  of  the  Michels  on  interferometer  to 
the  fringe  counter.  At  this  point,  it  splits  into  two  branches.  One  branch 
goes  to  the  sampler  where  the  signal  controls  the  samp’ 5  interval.  The 
second  branch  goes  to  the  sequencer  where  the  signal  ton tvole  the  frequency 
of  the  cosine  functions. 

The  summation  expressed  by  equation  (46)  it  performed  in  a  digital 

averager  or  summation  unit.  Signals  of  the  form  G(t  .  )  cosu>  t  are 

nn  n 

digitized  at  1024  ordinates,  and  each  of  these  values  is  stored  in  an  address 
of  the  memory.  When  the  mth  sample  point  is  taken  from  the  interferogram, 
the  term  G(t^^)  cosw^t  is  produced  by  the  multiplier.  This  is  digitized 
and  added  to 


nsm-1 

2  G(t  }  cob«  t 
,  nn  n 
nsl 

which  is  already  stored  in  the  memory.  The  previous  data  are  read  out  tx 
each  address,  modified  by  the  new  data,  and  the  result  is  rewritten  into 
the  address.  When  this  has  been  done  for  all  the  addresses,  the  Bununatlr  . 
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S  G{t  .  )  cosw  t 
n»l  nh  n 

is  stored  in  the  memory.  Corresponding  ordinates  of  the  successive 
signals  are  stored  in  the  same  address.  One  additional  signal  per  second 
may  be  added  to  the  memory.  An  external  trigger  initiates  the  digitization 
storage  cycle,  and  an  external  signal  advances  the  memory  from  one 
address  to  the  next. 

The  data  in  the  memory  are  monitored  on  an  oscilloscope  and  may  be 
plotted  on  an  X-Y  recorder  as  soon  as  the  interferogram  is  obtained. 

As  additional  terms  are  added  to  the  synthesized  spectrum,  greater 
detail  appears  and  the  resolution  increases.  Figures  37  and  38  show 
the  effect  of  additional  terms  on  the  synthesised  spectra. 

The  synthesized  spectra  repeat  at  intervals  of  1  /h.  The  external 
advance  is  adjusted  to  limit  the  repeating  spectra  to  the  last  24  addresses 
of  die  memory.  The  non-repeating  spectra  are  contained  in  the  first  1000 
address  points. 


45 


Description,  of  the  Individual  Circuits 

The  foregoing  section  discussed  the  operation  of  the  complete 
instrument  while  this  section  is  devoted  to  a  detailed  description  of  some 
of  the  more  important  components  and  circuits.  The  first  of  these  is  the 
zero  path  difference  detector  which  initiates  the  automatic  operation  of 
the  system. 

Zero  path  difference  detector 

It  has  been  shown  that  the  spectrum  may  be  obtained  by  synthesizing 
the  cosine  transform  of  the  interferogram.  This  requires  an  accurate 
knowledge  of  the  point  x  »  0  or  t  «  0,  depending  upon  which  function,  F(x) 
or  G(t)t  is  being  transformed.  Errors  in  the  location  of  this  point  perturb 
the  shape  of  the  scanning  function,  which  becomes  asymmetrical  and  has 
important  negative  minima*  Part  of  a  typical  interferogram  and  the  effect 
of  errors  in  the  determination  of  the  zero  point  are  shown  in  Figure  14  . 
Figure  15  shows  the  effect  of  a  small  error  on  the  shape  of  the  scanning 
function.  The  term  e  /x  is  a  measure  of  the  phase  error  where  *  is  the 
error  ana  ^  is  the  wavelength  of  the  radiation  [25].  U  the  error  is  small 
with  respect  to  the  shortest  wavelength,  the  scanning  function  changes 
only  slightly  as  a  function  of  wavelength. 

When  broad  spectra  are  studied,  the  interferogram  in  the  region  of 
zero  path  difference  contains  only  a  few  fringes  of  rapidly  decreasing 
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intensity*  The  location  of  the  point  of  zero  path  difference  can  be  deter¬ 
mined  with  precision.  An  electronic  circuit  could  select  this  point  by 
properly  operating  upon  the  interferogram  signal.  However*  for  narrow 
band  spectra*  the  region  near  aero  path  difference  has  many  fringes  of 
nearly  constant  intensity.  In  this  case*  an  uncertainty  of  several  fringes 
may  exist  and  an  electronic  circuit  is  not  able  to  select  the  proper  point. 
Examples  of  typical  interferograms  are  presented  in  Figure  16  . 

The  zero  path  difference  detector  consists  of  an  additional  optical 
system  and  the  associated  electronics.  The  optical  system  is  clamped  to 
the  sides  of  the  lamellar  grating.  A  lamp  and  collimator  produce  a 
narrow  beam  of  light  which  strikes  the  grating  at  near  grazing  incidence 
and  covers  a  narrow  section  of  all  the  facets.  Stops  isolate  the  zeroth 
order  interference  pattern*  and  a  lens  focuses  it  onto  the  detector.  The 
optical  system  and  the  mounting  details  are  shown  in  Figure  34  •  In  the 
figure*  K  indicates  the  light  source  and  collimator  while  L>  is  the 
detector. 

Figure  17  shows  the  configuration  of  the  U.  .ellar  grating  for  the 
general  case  and  for  the  case  of  zero  path  diLurence.  As  a  result  of 
constructive  interference  between  rays  1  and  2,  the  output  of  the  detector 
is  a  maximum  when  the  groove  depth  is  zero.  The  grating  then  behaves 
like  a  plane  mirror.  The  angle  of  incidence  of  the  light  beam  1b  critical 
and  must  be  carefully  adjusted  to  obtain  the  response  shown  in  Figure  18  • 
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This  figure  shows  the  detector  response  in  the  region  of  zero  path 
difference  as  the  movable  facets  are  driven  in  both  directions, 

A  block  diagram  and  schematic  of  the  electronic  system  is  shown  in 
Figure  19  •  The  triangle  with  an  "A"  is  used  to  designate  an  operational 
amplifier,  and  this  notation  will  be  used  whenever  operational  amplifiers 
are  employed. 

Input,  output  and  intermediate  waveforms  are  shown  in  Figure  18  , 
The  input  is  amplified  and  clipped  to  isolate  the  peak  which  is  then  differ¬ 
entiated  with  the  resulting  signal  passing  through  zero  when  the  input 
signal  reaches  the  peak.  An  amplifier  and  clipper  circuit  are  used  to 
shape  the  differentiated  input  signal  into  a  rectangular  pulse*  The  trailing 
edge  of  this  pulse  corresponds,  in  time,  to  the  peak  of  the  detector 
response,  A  narrow  voltage  spike  is  generated  by  differentiating  the 
rectangular  pulse.  The  spike  triggers  a  circuit  which  closes  the  relay. 
This  completes  the  signal  paths  used  by  the  outputs  of  the  two  interfer¬ 
ometers,  The  spike  may  be  advanced  or  retarded  slightly,  with  respect 
to  the  peak  of  the  detector  signal.  This  permits  a  fine  adjustment  of  the 
detected  point  of  zero  path  difference.  The  detected  point  of  zero  path 
difference  is  compared  to  the  actual  value  in  Figure  20  .  The  coincidence 
of  the  actual  and  detected  points  of  zero  path  difference  is  repeatable  if 
the  grating  is  driven  in  one  direction.  The  two  points  do  not  coincide 
when  the  grating  is  driven  in  the  opposite  direction.  This  is  not  a 
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limitation  to  the  system  performance  since  it  is  only  necessary  to  take 
data  from  x  =  0  to  x  =  Df  where  D  is  the  maximum  path  difference. 

Input  amplifier  and  sample  and  hold  circuit 

The  interferogr^pi  function  was  defined  in  equation  (16)  as 

CO 

F(x)  =  E(x)  -  yE{0)  a  f  E(v)  cos2irvxdv  (48) 

V,'’ 
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The  expression  jE(0)  is  the  asymptotM  value  approached  by  E(x)  as  the 
path  difference  becomes  large.  It  is  a  dc  value  which  must  be  removed 
or  an  error  will  be  produced  in  the  synthesized  spectrum.  A  network 
is  connected  to  the  input  amplifier  to  cancel  this  dc  value.  The 
cancellation  is  done  experimentally  by  adjusting  the  dc  level  of  the 
input  to  the  sample  and  hold  circuit  to  zero.  A  standard  analog  sample 
and  hold  circuit  is  used  in  the  system.  This  supplies  the  dc  input  to 
the  multiplier.  Figure  21  „ 


Multiplier 

The  system  uses  a  quarter  square  analog  multiplier.  This  type 
was  chosen  because  it  exhibits  superior  phase  and  frequency  response 
in  the  region  of  interest,  <2500  cps.  It  uses  biased  diode  squaring 
network^  and  these  are  the  cause  of  moderate  waveform  distortion  when 
the  dc  input  becomes  less  than  1  volt.  Compensation  is  accomplished 


by  means  of  a  level  detector  and  a  pair  of  amplifiers.  One  of  these 
increases  the  dc  level  by  a  factor  of  10  while  the  other  alternates  the 
output  of  the  multiplier  by  the  same  factor*  The  multiplier  then  operates 
over  a  more  linear  portion  of  its  characteristic.  An  absolute  value 
circuit  and  level  detector  switches  these  amplifiers  in  and  out  of  the 
signal  path  as  required*  The  multiplier  and  its  compensation  system 
are  shown  in  the  complete  system  diagram.  Figure  27  . 

Frequency  synthesizer  and  gate 

The  frequency  synthesizer*  gate*  and  gate  control  circuitry  are  used 
to  generate  cosine  functions  of  the  proper  frequency.  The  frequency 
synthesizer  is  a  standard  unit  with  provision  for  remotely  controlling  the 
frequency.  The  frequencies  are  synthesized  by  the  combination  of  various 
frequencies  which  are  derived  from  a  crystal  controlled  master  oscillator. 
Thus*  the  accuracy  of  these  frequencies  is  precisely  controlled.  The 
proper  frequency  is  selected  by  a  switching  network  consisting  of  four 
10-contact  rotary  step  switches.  The°“  are  driven  by  the  sequencer. 
Frequencies  may  be  selected  in  1  cycle  steps  from  1  ops  to  9999  cps. 

The  frequency  response  of  the  synthesizer  is  flat  over  the  range  of  intere  ■ , 
<2,5  kc,  and  varies  only  0*5%  from  0  to  10  kc.  Harmonics  are  down  at 
least  40  db. 

A  standard  type  of  gate  circuit  is  used*  it  is  operated  by  a  pulse 
generator  which  is  triggered  by  the  output  of  the  90°  detector. 


90- degree  detector 


Cosine  functions  are  obtained  from  the  frequency  synthesizer  by- 
means  of  a  90-degree  detector  and  a  gate  circuit#  The  90°  detector 
opens  the  gate  at  the  peak  of  the  sinusoidal  signals.  This  permits  a 
cosine  function  to  pass  through  the  gate#  The  gate  automatically  closes 
?fier  a  one -second  time  delay.  The  90-degree  detector  should  meet  the 
following  specifications: 

1.  It  should  be  able  to  detect  the  90°  point  of  a  sine  wave  while 
operating  over  a  frequency  range  of  more  than  3  decades,  from  1  cps  to 
2500  cps. 

2.  It  should  not  require  tuning  or  readjustment  as  the  frequency  is 
varied. 

3.  The  output  of  the  frequency  synthesizer  should  be  the  only  input 
signal  required. 

Many  standard  techniques  are  available  for  triggering  circuits  at 
perpendicular  points  on  a  wave  form.  These  generally  use  a  level 
detector  and  some  phase  shifting  scheme.  None  of  these  will  operate 
satisfactorily  for  this  particular  application.  Level  detectors  themselver 
are  unstable  in  the  region  of  the  peak  of  a  sine  wave.  Phase  shifters  used 
in  conjunction  with  level  detectors  do  not  have  the  desired  bandwidth. 
Phase  splitters,  such  as  are  used  in  single  sideband  equipment,  have 
neither  the  required  bandwidth  nor  the  proper  amplitude  response.  The 
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heterodyne  and  phase  shifting  technique  ha#  also  proven  unsatisfactory 
for  this  application.  Shifting  phase  by  differentiation  or  integration  is 
not  practical  when  the  bandwidth  extends  over  three  decades. 

The  technique  that  was  developed  is  to  utilize  the  sinusoidal  output 
of  the  frequency  synthesizer  and  shape  it  into  a  series  of  very  narrow 
pulses.  It  is  important  to  maintain  coincidence  between  the  peak  of  the 
input  and  the  peak  of  the  pulse.  The  pulse  shaping  is  done  with  non-linear 
circuits.  Each  stage  generates  a  narrower  pulse.  The  last  stage  of  the 
shaping  circuitry  drives  a  level  detector.  Its  output  is  now  stable  in  the 
region  of  the  peak  of  the  input,  since  it  is  operating  on  a  narrow  pulse 
which  has  a  rapid  rise  and  fall.  The  wave-shaping  circuits  are  the  key 
to  the  successful  operation  of  the  90-degree  detector. 

A  diagram  of  the  shaping  circuit  is  presented  in  Figure  22  ,  A 
schematic  of  two  of  the  shaping  networks  and  their  amplifiers  is  shown  in 
Figure  23  .  The  output  of  the  frequency  synthesizer  is  amplified,  and 
the  resulting  signal  goes  to  the  gate  circuit  and  also  to  the  first  stage  of 
the  shaping  circuit.  This  is  a  non-linear,  biased-diode  network.  The 
component  values  have  been  chosen  to  produce  an  approximately  quartic 
characteristic.  Thus,  the  output  current  is  related  to  the  input  voltage 
by 
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where  k  is  a  constant  chosen  to  give  the  components  reasonable  values. 

Six  line  segments  are  used  in  the  approximation. 

The  output  of  the  first  shaping  network  drives  the  second  one  directly. 

5 

This  has  an  approximate  quintic  characteristic,  iQ  »  k  E  ,  as  have  all  the 
remaining  networks. 

With  reference  to  Figure  23  ,  the  output  current  of  the  biased  diode 
network  flows  through  the  2  K  resistor  and  develops  a  voltage  which  is 
then  amplified.  The  amplifier  also  serves  tc*o  additional  functions.  It  is 
a  buffer  between  the  networks  to  prevent  loading  effects,  and  it  permits 
adjustment  of  the  dc  level.  The  ability  to  adjust  the  dc  level  is 
important  to  the  proper  operation  of  the  circuit,  and  each  stage  is  provided 
with  a  network  by  which  a  dc  voltage  may  be  added  to  or  subtracted 
from  the  signal. 

The  final  shaping  circuit  is  a  non-linear  amplifier.  A  field  effect 
transistor  is  used  in  the  feedback  loop  of  an  operational  amplifier  to 
provide  the  effect  of  a  voltage  variable  resistor.  The  connection  of  the 
field  effect  transistor  is  shown  in  Figure  22  ,  In  this  case,  the  effective 
resistance  of  the  transistor  increases  as  the  output  voltage  increases. 

Since  the  transfer  function  of  the  amplifier  is 

6o  ^feedback  . 
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in  input 


the  gain  increases  aonlinsarly.  The  net  effect  of  this  circuit  is  to  produce 
very  narrow  pulses  which  occur  at  the  peaks  of  the  original  input  waveform. 
The  output  of  the  nonlinear  amplifier  drives  a  level  detector,  which  is 
now  stable  in  the  region  of  the  peak  of  the  input  sinusoid,  since  it  operates 
on  a  narrow  pulse. 

Figure  24  shows  the  waveforms  obtained  at  various  points  in  the 
circuit.  It  is  observed  that  each  succeeding  stage  generates  a  narrower 
pulse.  The  output  of  the  level  detector  is  shown  for  two  values  of  the 
sweep  speed.  In  this  case,  the  level  detector  is  triggered  1.44°  before 
the  90°  point  is  reached.  The  original  circuit  design  did  not  include  the 
nonlinear  amplifier  and,  in  this  case,  the  stable  triggering  point  was  6° 
before  the  peak.  Since  the  pulse  is  very  narrow  but  has  a  finite  width, 
the  triggering  point  will  be  slightly  ahead  of  the  90°  point  unless  phase 
shift  elsewhere  in  the  circuitry  causes  a  delay. 

The  accuracy  of  the  circuit  is  dependent  upon  the  initial  adjustment, 
the  amplitude  stability  of  the  various  signals,  phase  shift,  and  power 
supply  fluctuations.  Gain  controls  and  a  dc  balancing  network  are 
provided  for  each  stage.  The  initial  adjustment  of  these  controls  is 
critical  because  of  the  nonlinear  circuits,  high  amplification,  and  the 
saturation  characteristics  of  the  amplifiers.  Once  the  circuit  has  been 
properly  adjusted,  no  further  adjustments  are  required  as  the  frequency 
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The  biasing  networks  employ  zener  diodes  to  stabilize  the  operation 
of  the  shaping  networks,  A  stable  power  supply  with  low  ripple  voltage 
is  very  important  to  the  proper  operation  of  the  circuits* 

A  combination  of  phase  shift  in  the  shaping  networks  and  lack  of 
sufficient  high  frequency  response  in  the  amplifiers  causes  the  operation 
of  the  circuit  to  deteriorate  beyond  2500  cps.  This  range  of  operation  is 
adequate  for  the  present  system. 

Figure  25  shov/e  a  plot  of  the  measured  triggering  points  for  the 
level  detector  as  a  function  of  frequency.  Triggering  error  is  less  than 
2  degrees  in  the  range  from  2  cps  to  3  ke.  This  is  a  considerable 
improvement  over  several  other  methods*  Level  detectors  become 
unstable  at  ubout  15  degress  from  the  peak  of  a  sinusoid.  Single  sideband 
pha^e  splitters  generate  2  signals  90°  HH  1. 5°  out  of  phase  with  each  other 
over  a  bandwidth  of  1  decade,  300-3000  cps,  and  the  amplitude  response 
is  not  flat.  Figure  25  also  shows  the  phase  shift  at  the  output  of  ths 
fourth  amplifier  as  compared  to  the  input,  as  a  function  of  frequency. 

Summation  unit 

A  standard  digital  summation  unit  is  used  to  sum  the  G(t^)  cosu^t 
terms  as  they  are  produced  in  the  multiplier.  The  output  of  the  multiplier 
is  digitized  at  1024  ordinates,  and  sach  of  these  values  is  stored  in  a 
separate  address  of  ths  memory.  Two  external  control  signals  are  used, 
Ths  first  is  derived  from  the  gating  circuit  which  produces  a  pulse  at  tbA 
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instant  the  gate  opens*  This  pulse  initiates  the  digitisation  and  sweep  cycle 
in  the  memory.  The  second  control  signal  is  derived  from  a  fixed 
frequency  output  of  the  frequency  synthesiser.  It  ie  used  to  externally 
advance  the  memory  from  one  address  to  the  next*  These  two  control 
signals  ensure  that  the  digitiser  begins  operation  at  the  same  point  (the 
peak)  of  each  succeeding  cosine  function*  and  that  data  for  corresponding 
ordinates  are  stored  in  the  same  address.  The  synthesised  spectrum  will 
repeat  at  intervals  of  1/h.  The  highest  frequency  in  the  synthesized 
spectrum  is  1  /Zh,  and  the  repetition  rate  of  the  external  advance  signal 
is  adjusted  to  scale  this  frequency  to  the  1000  address.  The  last  24 
addresses  then  carry  part  of  the  repeated  spectrum,  ft  is  important  to 
properly  choose  the  sampling  intervals  to  reduce  the  effect  of  overlap  of 
noise  frequencies.  The  techniques  for  this  are  discussed  in  the  first 
section  of  this  report. 

Fringe  counter 

The  sinusoidal  reference  signal  is  derived  from  the  Micheleon 
interferometer  and  is  used  to  control  the  sampling  interval  and  the 
frequency  of  the  cosine  functions.  The  zero  crossings  are  used  as  the 
reference  points*  and  each  cycle  indicates  a  0.27305  micron  additional 
displacement  of  the  movable  facets.  The  sampling  theorems  do  not 
require  that  samples  be  taken  this  closely*  if  the  shortest  wavelength  is 


10  microns.  Therefore,  a  counter  is  used  which  enables  samples  to  be 
taken  at  each,  every  other,  every  or  16^  cycle. 

The  sinusoidal  reference  signal  is  converted  to  pulses  and  delayed 
in  a  variable  delay  circuit.  The  delay  circuit  consists  of  a  single-shot 
pulse  generator  with  a  variable  pulse  width.  The  delay  is  incorporated 
to  ensure  that  beginning  of  the  sampling  interval  and  the  output  of  the 
xero  path  detector  coincide.  As  a  result  of  mechanical  imperfections  in 
the  system,  and  vibration,  the  coincidence  cannot  be  repeated  more  than 
a  few  times  before  an  adjustment  is  required.  The  maximum  error 
introduced  if  the  delay  is  not  used,  ox  if  it  is  misad justed,  is  half  a 
mercury  fringe  or  •  136525  microns.  If  the  infrared  radiation  is  10 
microns,  the  phase  error  is  «  /x  *  .0136525.  With  reference  to  Figure  16 
the  scanning  function  is  only  slightly  altered  and  no  appreciable  error 
results. 


The  fringe  counter  provides  the  "command"  signal  for  the  sample 
and  hold  circuit  and  the  sequencer.  The  sequencer  consists  of  counters, 
pulse  generators,  transistor  driving  circuits  and  four  rotary  step  switches. 
A  block  diagram  of  the  sequencer  is  shown  in  the  lower  right-hand  part 
of  Figures  ?.l  and  27  •  The  block  marked  "SSMV"  contains  the  pulse 
generator.  The  leads  which  control  the  frequency  generated  by  the 
frequency  synthesiser  are  connected  to  the  10  contact  Step  switches. 
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Each  time  a  pulse  enters  the  sequencer,  the  "units"  step  switch 
advances  one  contact  and  the  frequency  of  the  synthesiser  advances  one 
cycle  per  sec  end.  The  remaining  switches  operate  in  the  proper  counting 
sequence,  and  the  frequencies  can  be  advanced  from  1  epa  to  9999  cps  in 
steps  of  one  cps.  All  the  step  switches  may  be  reset  to  aero  by  means  of 
a  "homing"  switch* 

The  coils  of  the  step  switches  are  energized  by  the  application  of 
a  pulse  from  the  transistor  driving  circuits*  No  switching  takes  place, 
however,  until  the  coils  are  de-energized*  This  introduces  a  delay  of 
approximately  40  milliseconds  between  the  application  of  the  switching 
pulse  and  a  change  in  the  frequency  of  the  frequency  synthesizer. 

A  70  millisecond  delay  is  introduced  into  the  gate  control  circuitry 
to  ensure  that  the  gate  epens  after  the  proper  frequency  has  been 
established  and  the  switching  transients  have  died  out* 

Control  circuitry 

Figure  26  presents  a  timing  diagram  which  shows  the  important 
relationships  between  the  several  inputs,  outputs  and  control  signals*  A 
complete  block  diagram  of  the  electronic  system  is  shown  in  Figure  27  • 
This  shows  the  control  circuits  in  detail* 

The  counter  signal,  comes  from  the  fringe  counter  and  is  derived 
from  the  reference  input.  Each  time  a  negative  transition  occurs,  the 
sampler  takes  a  new  sample  and  the  output  of  the  frequency  synthesizer 
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i*  advanced  one  cycle  per  second.  The  delay,  t  u  the  40-millisecond 
switching  delay  between  the  application  of  the  switching  pulse  and  the 
actual  switch  closure.  The  90-degree  detector  generates  a  pulse  each 
time  the  output  of  the  frequency  synthesizer  passes  through  a  maximum. 
The  signal  which  controls  the  sampling  in  the  sample  and  hold  circuit  is 
also  available  aa  an  output  of  that  circuit.  This  signal  is  delayed  70 
milliseconds  and  fed  to  one  input  of  a  bistable  circuit.  The  signal  from 
the  90-degree  detector  is  fed  to  the  other  input.  The  sample  and  hold 
activation  signal  sets  the  output  transistor  to  the  "off"  state.  The 
90-degree  detector  resets  it  to  the  "on"  state.  This  creates  a  negative 
voltage  transition  which  operates  the  gate-driver  circuit  and  opens  the 
gate.  It  also  produces  a  trigger  signal  which  begins  the  digitizing  sequence 
in  the  summation  unit.  The  gate  is  closed  after  a  predetermined  time 
delay  by  the  gate  driving  circuit.  When  the  output  of  the  gate  enters  the 
multiplier,  it  is  multiplied  with  the  value  of  the  interfcrogram  sample 
and  is  immediately  digitized  and  stored  in  the  summation  unit. 


TEST  RESULTS  AND  THE  SYNTHESIZED  SPECTRA 
OBTAINED  WITH  THE  INSTRUMENT 


Experimental  Procedure 

The  completed  real-time  Fourier  transform  synthesiser  was 
subjected  to  a  series  of  tests  designed  to  evaluate  its  performance* 

Since  it  will  be  used  to  synthesize  the  spectra  produced  by  both  broad 
and  narrow  spectral  band  sources,  known  functions  with  these  types  of 
spectra  were  used.  The  entire  system  was  not  operated  during  those 
initial  tests.  The  zero  path  difference  detector  was  disconnected,  and 
the  relay  it  controls  was  manually  operated. 

The  performance  in  the  synthesis  of  continuous  spectra  was 
investigated  by  utilizing  single  pulses  of  various  widths.  Rectangular, 
triangular  and  cosine  pulses  were  used.  They  were  derived  from  a  low- 
frequency  signal  generator,  and  an  identical  generator  supplied  die 
reference  signal.  These  signals  were  inserted  into  the  normal  data  inpu 
and  their  amplitudes  were  adjusted  to  simulate  the  normal  input  levels. 
The  signals  were  monitored  on  a  strip  chart  recorder,  and  the  sampling 
sequence  was  manually  initiated  at  the  proper  instant.  The  summatio;' 
of  £(t^)  cosw^t  terms  was  stopped  at  the  proper  time  by  switching  the 
summation  unit  from  die  "count"  mode  to  the  "display"  mode. 

The  test  functions  will  be  discussed  mathematically  using  Bracewa;  i 
[26]  notation.  This  notation  is  now  defined. 


7  * 

The  rectangular  function  of  unit  height,  area,  and  base  is  defined 
as  II(x)« 


H(x) 


li 


0  for  ix.  >  i 
1  for  x;  <| 
J  for  fx,  a  | 


The  triangular  function  of  unit  height  and  area  is  defined  as  (x)* 


(*) 


^  0  for  |  x  i  >  1 
|^1  -  ix}  for  IX!  <  1 


This  function  is  also  the  t.  elf -convolution  of  H(x). 


The  sampling  or  replicating  symbol  IIX(x)  is  defined  as 


in(x)  =  2  A(x-n) 
n=-cc 


T 
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The  rectangular  functions  of  unit  height  may  be  represented 
mathematically  as  H(t/b),  where  t  is  time  and  b  is  width.  In  a  similar 
manner,  the  triangular  functions  become  (t/b/2)  and  the  cosine  pulses 
are  expressed  as  I£{t/b)  cosir(t/b).  Figure  28  shows  the  rectangular 
functions  and  both  their  synthesized  and  theoretical  spectra.  Figures  29 
and  30  show  similar  data  for  the  triangular  and  cosine  pulses, 
respectively.  The  fluctuations  in  the  synthesized  data  are  due  to 
digitization  error  and  recorder  jitter.  Figure  31  shows  a  set  of 
normalized  error  curves  for  the  preceding  data.  The  error  is  normal;ze.i 
to  the  amplitude  of  the  dc  component  of  the  spectrum,  and  its  absolu'  a 
value  is  plotted.  The  greatest  error  occurs  for  shorter  pulses.  This 
is  so  for  .wo  reasons:  First,  only  a  few  terms  make  up  the  synthesized 
spectrum  and  an  insufficient  number  are  used  to  obtain  good  convergence 
to  the  true  spectrum.  Second,  because  of  the  sampling  rate  limitations 
imposed  by  the  components  of  the  system,  the  fastest  allowable  sampling 
rate  is  insufficient  to  allow  sampling  at  twice  the  highest  frequency  of  the 
spectrum.  The  spectra  of  these  signals  theoretically  extends  to  infinity, 
but  amplitudes  soon  become  negligible  for  all  practical  purposes.  The 
region  of  the  spectrum  in  which  the  amplitudes  are  non-negllgible  is 
wider  for  short  pulses  than  it  is  for  long  ones.  Thus,  the  allowable 
sampling  rate  is  not  fast  enough  to  prevent  overlap  of  spectral  componer>' 
with  non -negligible  amplitudes. 
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The  accuracy  improves  for  long  pulses  because  a  greater  number 
of  samples  produces  a  larger  number  of  terms.  The  spectrum  estimate, 
therefore,  more  closely  approaches  the  true  spectra.  Also,  the  non" 
negligible  portion  of  the  spectrum  is  narrower  and  the  overlap  of  repeated 
spectra  occurs  for  components  with  much  smaller  amplitudes. 

The  spectrum  for  a  short  rectangular  pulse  was  synthesised  from 
only  three  samples  while  the  short  cosine  and  triangular  pulses  were 
synthesized  from  only  four  samples.  Sampling  rate  limitations  are 
discussed  further  in  the  next  section. 

The  preceding  functions  have  continuous  spectra  which  contain 
no  rapid  variations  or  discontinuities.  It  is  therefore  possible  to  syn¬ 
thesis  e  the  spectra  with  a  finite  number  of  terms  and  obtain  good  accuracy, 
provided  the  number  of  terms  is  sufficient*  This  is  not  the  case  when  the 
spectrum  is  discontinuous,  since  Gibbs  phenomenon  occurs.  An  infinite 
number  of  terms  produces  a  synthesized  spectrum  which  exhibits  "ringing" 
at  the  dis continuities. 

Let  us  now  investigate  the  synthesized  spectrum  estimate  when 
functions  having  discontinuous  spectra  are  involved.  For  this  purpose, 
finite  trains  of  rectangular  and  triangular  pulses,  and  a  truncated  cosine 
wave,  were  used. 

Mathematically,  the  rectangular  pulse  train  can  be  represented  as 


f(t)  a  n(t/b){{III(t/e)*2lI<t/a)}  -  1} 


where  *  denotes  convolution  and  the  parameters  are  defined  by  the 
following  diagram. 


The  rectangular  pulse,  convolved  with  the  replicating  function  produces 
a  pulse  train.  The  dc  value  is  removed  and  the  result  is  truncated  by 
a  second  rectangular  function.  In  the  frequency  domain,  the  spectrum  is 
the  spectrum  of  the  rectangular  pulse,  multiplied  by  the  sampling  function, 
which  is  then  convolved  with  the  scanning  function.  In  this  case,  the 
scanning  function  is  the  Fourier  transform  of  the  truncating  function. 


Thus,  the  line  spectrum  of  an  infinite  pulse  train  is  convolved  with  a 

scanning  function  of  the  form  -81^  irt  =  sine  t.  The  synthesized  spectrum 

irt 

will  be  a  series  of  sine  functions  located  at  the  harmonics  of  the  fundamen¬ 


tal  frequency. 

The  truncated  cosine  wave  is  represented  by 


f(t)  =  II(t/b)  cos  wt 


G2 

The  spectrum  is  simply  a  pair  of  impulse  functions  convolved  with  the 
scanning  function.  In  this  case,  one  of  the  impulse  functions  occurs  for 
negative  frequencies  which  have  no  meaning. 

The  triangular  function  may  be  represented  in  the  general  case  as 


f<t)  =  U{t/b)[{m<t/c)*  2  (t/a) )  -  1] 

b - 
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Figure  32  shows  the  waveforms  and  the  synthesized  spectra.  Each 
waveform  was  sampled  at  2-second  intervals,  and  500  samples  were 
obtained.  Table  1  shows  the  measured  and  theoretical  values  for  the 
spectral  components. 

Frequencies  may  be  determined  with  an  error  of  less  than  2%.  The 
accurate  determination  of  amplitudes  is  somewhat  more  difficult.  This 
is  due  to  the  fact  that  there  are  only  a  finite  number  of  addresses  in  which 
data  are  stored.  The  resolution  of  the  system  improves,  and  the  peak  of 
the  synthesized  harmonic  may  actually  fall  between  two  address  points 
when  many  samples  are  taken.  In  this  case,  the  true  peak  will  not  be 
recorded.  This  is  true  of  any  system  which  utilizes  a  finite  number  of 
input  or  output  points.  The  oscilloscope  display  was  expanded,  and  the 
main  prak  was  found  to  be  stored  in  only  four  addresses.  The  peak  was 
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Table  1.  Result*  of  the  synthesis  of  the  spectra  of  truncated 
waveforms.  * 


Cosine  Wave 

Harmonic 

1 

3 

5 

Amplitude 

Theoretical 

1 

0 

0 

Synthesized 

1 

0 

0 

|  Error  j 

0 

0 

0 

F  requency 

Theoretical 

1 

0 

0 

Synthesized 

1 

0 

0 

(Error  | 

0 

0 

0 

Triangular  Wave 

Harmonic 

1 

3 

5 

Amplitude 

Theoretical 

1.0 

0.111 

.04000 

Synthesized 

1.0 

0. 1 094 

. 0i 885 

'Error  j 

0 

1.45% 

52.  8% 

Frequency 

Theoretical 

1.0 

3.00 

5.00 

Synthesized 

1.0 

3.05 

5.  10 

’Error  ’ 

0 

1.6% 

2.0% 

? 


0 

0 

0 


0 

0 

0 


7 


.  021 


9 

0 

0 

0 


0 

0 

0 


9 


.01235 


7.  00 


9.  00 
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Table  l.  (Continued), 


Rectangular  Wave 

Harmonic 

1 

3 

5 

7 

9 

Amplitude 

Theoretical 

1 

-0.  333 

0.200 

-C.  1428 

0.  Ill 

Synthesized 

1 

-0.391 

0.221 

-0. 1465 

0.129 

|  Error j 

0 

18.  5% 

10.  5% 

2.6% 

16.4% 

Frequency 

Theoretical 

1 

3.00 

5.00 

7.  00 

9.00 

Synthesized 

1 

3.00 

5.00 

7.10 

9.04 

Error 

0 

0 

0 

1.428% 

0.445% 

All  the  data  are  normalized  to  the  amplitude  and  frequency  of  the 
fundamental. 
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not  stored.  This  causes  an  amplitude  error  in  the  synthesized  harmonic. 
Methods  of  reducing  this  error  are  discussed  in  the  following  section. 

A  test  of  the  resolution  of  the  system  was  made  with  these  signals. 

The  theoretical  resolution  is  .  0005  cps  while  the  measured  resolution  is 
approximately  the  same.  Another  rest  was  made  with  a  triangular  wave¬ 
form  and  61  samples  were  taken.  The  measured  and  theoretical  resolutions 
were  nearly  the  same  again,  but  the  amplitudes  of  the  harmonic  components 
were  more  nearly  equal  to  the  theoretical  values.  The  peaks,  representing 
harmonics,  were  stored  in  more  addresses,  and  the  actual  location  of  the 
peak  was  more  closely  determined. 

The  real-time  Fourier  transform  synthesizer  was  evaluated  under 
actual  operating  conditions.  The  experimental  apparatus  is  shown  in 
Figure  33  .  Figure  34  is  a  photograph  of  the  optical  system  while 
Figure  35  is  a  photograph  of  the  equipment.  In  Figure  34  ,  E  is  the 
lamellar  grating,  K  and  L  are  the  optics  for  the  zero  path  difference 
detector,  G  is  the  globar,  P  is  the  monochromator,  and  PM  is  the  parabolic 
mirror. 

The  lamellar  grating  was  illuminated  with  parallel  radiation  produced 
by  the  globar  which  is  at  the  focal  point  of  the  parabolic  mirror.  After 
reflection  from  the  face  of  the  grating,  this  radiation  is  focused  by  the 
optical  system  into  the  monochromator.  This  was  used  as  both  an  optical 
filter  and  a  detector.  The  output  of  the  monochromator  was  amplified  and 


EXPERIMENTAL  SETUP 


OPTICAL  SYSTEM 
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LAMELLAR  GRATING  L  =  DETECTOR 

GLOBAR  P  =  MONOCHROMATOR 

LIGHT  SOURCE  AND  COLLIMATOR  PM  =  PARABOLIC  MIRROR 
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EXPERIMENTAL  SET-UP 


FIGURE  35 
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connected  to  the  it^ut  of  the  Fourier  transform  synthesizer.  The  inter- 

ferogra-  i  recorded  on  a  strip  chart  recorder.  The  reference  signal 

was  obtained  from  the  Michels  on  interferometer,  and  the  zero  path 

difference  detector  supplied  the  initiating  pulse. 

Figure  36  shows  the  interferogram  and  the  corresponding  synthesized 

spectrum.  The  optical  path  difference  was  varied  from  -400  to  +400 

microns.  The  samples  were  taken  over  the  0  to  400  micron  region.  It 

required  approximately  200  seconds  to  ob*a*n  the  data.  The  synthesize 

spectrum  is  always  monitored,  and  Figure  3 7  shows  the  increase  in 

detail  as  the  path  difference  is  increased. 

The  interferogram  is  produced  by  a  narrow  band  of  radiation  with 

a  central  wavenumber  of  approximately  605  cm”1.  This  corresponds  to 

a  wavelength  of  1/v^,  or  approximately  16.5  microns.  The  halfwidth  of 

the  synthesized  spectral  line  is  Av  =  30.5  cm  1  and  the  theoretical 

JLi 

resolution  of  the  instrument  is  A =  25  cm  1  for  the  path  difference 
utilized. 

Figure  38  shows  the  result  of  synthesizing  the  spectrum  produced 
by  the  interference  of  two  electrical  signals  of  slightly  different 
frequencies.  This  also  shows  the  increase  in  detail  and  resolution  as  thr> 
number  of  samples  (and  hence  terms)  increases. 


INTERFEROGRAM  AND  THE  SYNTH 
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Figure  37 
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Summary  of  the  data 

The  real-time  Fourier  synthesizer  has  been  tested,  under  actual 
operating  conditions.  It  is  able  to  synthesise  continuous  and  discrete 
spectra.  It  is  limited  in  the  maximum  frequency  at  which  samples  can 
be  obtained.  Thus,  greater  error  occurs  when  it  is  called  upon  to 
synthesize  short  pulses  with  broad  spectra.  Accuracy  improves  as 
the  pulses  become  wider  and  the  spectra  narrower.  The  theoretical  and 
measured  resolutions  were  approximately  the  same.  Accuracy  was 
satisfactory  when  it  was  used  to  synthesize  a  spectrum  from  the  corres¬ 
ponding  interferogram. 
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A  DISCUSSION  OF  ERRORS  AND  RECOMMENDED 
SYSTEM  IMPROVEMENTS 

The  synthesized  estimate  of  the  spectrum  will  deviate  slightly  from 
the  true  spectrum  because  of  the  following: 

1.  Errors  in  obtaining  the  interferogr&ms  are  unavoidable. 

2.  Errors  are  introduced  in  the  insertion  of  the  data  into  the 
computer. 

3.  Errors  arise  in  the  technique  of  computing  the  spectrum. 

4.  Errors  arise  due  to  limitations  of  the  component  parts  of  the 
Fourier  synthesizer. 

It  is  impossible  to  allow  the  path  difference,  when  taking  the  inter  - 
ferogram,  to  vary  from  to  +».  The  maximum  path  difference  is 
determined  by  the  limitations  of  the  interferometer.  The  interferogram 
is  therefore  truncated,  which  causes  an  unavoidable  error  and  limits  the 
resolution.  The  spectrum  estimate  is,  therefore,  the  true  spectrum 
convolved  with  the  scanning  function  of  the  instrument. 

Let  us  now  investigate  the  error  in  the  spectrum  when  an  improper 
interferogram  function  is  vised.  We  will  assume  a  noise-free  interfero¬ 
gram.  The  interferogram  function  has  been  defined  as 


F{x)  a  E(x)  =  fE(O) 
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where  E(x)  is  the  detector  response  and  f£(G)  is  its  asymptotic  value 
as  the  path  difference  is  made  large.  An  error  in  F{x)  is  made  if  jE{r>) 
is  improperly  removed  from  E(x).  For  an  error  of  AE,  the  error  in  the 
spectrum  will  be 


*(*)*£  = 4D^  <8i 

where  D  is  the  maximum  groove  depth  of  the  grating.  This  is  the  same 

as  the  Fourier  transform  of  a  rectangular  pulse  of  amplitude  AE  and 

width  4D.  The  error  adds  a  decaying  oscillation  with  a  period  in  v, 

to  the  spectrum.  This  error  is  ignored  in  several  digital  computation 

schemes  since  it  only  occurs  near  s  0  [22],  It  cannot  be  ignored  hero, 

for  all  wavenumbers  are  present  from  v  =  0  to  v  =  vw  where  v. ,  is  the 

M  M 

maximum  wavenumber.  Since  the  error  cannot  be  ignored,  it  must  bo 
removed,  The  input  amplifier  contains  a  balancing  network  to  remove 
|E(0)  from  E(x).  The  adjustment  is  made  experimentally  by  taking  the 
same  interferogram  several  times  and  adjusting  the  Ac  level  at  the 
input  of  the  sample  and  hold  circuit  to  zero.  Several  sets  of  data  are 
usually  required  to  reduce  AE  to  zero.  The  synthesized  spectrum  may  be 
monitored,  and  it  is  not  usually  necessary  to  obtain  the  entire  interfero¬ 
gram.  The  error  begins  to  appear  after  only  a  few  terms  in  the  synthesis 
have  been  obtained.  This  fact  saves  considerable  time  in  making  the 
adjustment. 
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The  spectrum  estimate  converges  to  the  true  spectrum  ae  the 
number  of  terms  in  the  Fourier  series  expansion  becomes  very  large, 
provided  the  spectrum  is  continuous.  A  finite  series  expansion,  produced 
by  an  interferogram  of  finite  length,  shows  spurious  oscillations  in  the 
spectrum  in  the  region  of  a  discontinuity.  The  amplitude  of  these 
oscillations  remains  constant  a  a  the  number  of  terms  increases,  but  the 
spurious  peaks  crowd  toward  the  discontinuity  [  23],  The  peak  separation 
is  given  by 


This  is  the  Gibbs  phenomenon. 

As  discussed  in  section  1,  Gibbs  phenomenon  may  be  removed  by 
a  process  known  as  "apodization,  "  in  which  the  interferogram  is  truncat' d 
hy  some  non-rectangular  function.  The  function  is  chosen  so  that  its 
Fourier  transform  (the  scanning  function)  has  a  main  lobe  and  low  side 
lobes.  A  common  apodizing  function  is  the  triangle  shown  in  Figure  5  . 

This  produces  a  scanning  function  of  the  form  (  -  —  )  while  rectangular 
truncation  produces  a  )  scanning  function.  The  (~~^)  function  has 

a  half  width,  for  the  same  number  of  terms,  which  is  twice  as  wide  as 

.  Thus,  for  triangular  apodization,  the  resolving  power  is  decrease^1 

I 

by  a  factor  of  two.  However,  the  scanning  function  for  triangular 
apodization  has  adjacent  side  lobes  which  have  amplitudes  about  4.  5%  of 


the  main  peak  and  they  do  not  oscillate.  The  expression  — ■ *  hac 
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adjacent  maximum  with  amplitudes  more  than  10%  of  the  main  peak.  It 
also  has  negative  valued  and  the  first  minimum  has  an  absolute  value  of 
approximately  20%  of  the  main  peak.  Other  apodizing  functions  can  be 
used,  and  the  choice  depends  on  the  characteristics  of  the  spectrum.  The 
present  system  does  not  contain  an  apodizing  system  since  an  additional 
multiplier  is  not  available.  This  feature  will  be  added  to  the  system  as 
soon  as  one  is  obtained. 

The  accuracy  of  the  real-time  Fourier  transform  synthesizer  may 
be  improved  by  changing  several  of  its  present  components.  Provision 
has  been  made,  in  the  design  of  the  system,  to  install  a  more  accurate 
analog  multiplier.  The  static  accuracy  of  the  present  unit  is  1%  while 
the  new  one  will  have  a  specified  accuracy  of  0. 2%.  The  installation  of  fie 
new  unit  will  free  the  present  one  so  it  may  be  used  in  the  apodization 
system. 

The  summation  unit  contains  an  analog  to  digital  converter.  It  is 
an  F*  system,  and  its  linearity  is  specified  as  .  5%.  It  digitizes  the 
average  value  of  the  signal  during  the  digitizing  interval.  This  causes 
the  frequency  response  to  begin  decreasing  for  frequencies  in  excess  of 
200  cps.  This  problem  could  be  solved  with  the  addition  of  a  separate, 
high  speed  analog  to  digital  converter,  with  provision  of  holding  the  si'jnr' 
value  during  the  digitizing  interval.  A  high  speed  unit  would  require  lers 
digitising  time,  and  these  facts  would  improve  the  frequency  response  of 
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the  converter -summation  unit  combination.  The  accuracy  of  digitization 
could  be  improved  by  using  a  converter  which  can  digitize  10  or  12  bits. 

This  gives  an  accuracy  of  0,  1%  or  0.  025%  respectively.  The  summation 
unit  can  accept  a  pulse  input  so  the  output  of  the  additional  converter  can 
be  fed  into  the  summation  unit  with  a  minimum  of  interfacing. 

In  section  3,  it  was  mentioned  that  the  point  where  the  peak  occurs 
in  the  discrete  line  spectra  may  not  be  recorded  due  to  a  finite  number  of 
channels.  This  is  a  limit  to  the  accuracy  of  resolving  amplitude  differ¬ 
ences.  This  is  present  in  any  system  in  which  only  a  finite  number  of 
points  are  recorded,  in  this  case  1024.  The  number  of  recorded  points 
can  be  increased  by  the  addition  of  a  larger  memory  with  2048  or  4096  <?■.’  Ire 
addresses.  A  solution  to  this  problem  involved  both  the  addition  of  more 
memory  addresses  and  an  interpolating  scheme  to  obtain  information 
located  between  address  points. 

The  maximum  allowable  sampling  rate  ie  limited  by  several  factors. 
First,  there  is  an  unavoidable  delay  in  switching  the  frequency  syntherizer 
with  stepping  switches.  Second,  it  is  required  that  at  least  a  half  cycle 
of  the  lowest  frequency  generated  by  the  frequency  synthesizer  be  stored 
in  the  memory.  The  lowest  frequency  is  presently  one  cycle  per  second.. 

The  maximum  sampling  rate  can  be  increased  about  five  times.  This 
requires  a  new  switching  scheme  and  increased  frequency  response 
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capabilities  of  the  memory  unit.  Thin  feature  will  be  added  as  soon  as 
an  additional  A-D  converter  become  available.  The  counting  circuits 
for  this  scheme  are  already  included  in  the  sequencer. 
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CONCLUSION 


Various  techniques  for  obtaining  the  Fourier  transform  have  been 
studied.  A  technique  was  developed  to  determine  the  transform  in  real 
time.  A  Fourier  transform  synthesiser  has  been  designed,  constructed 
and  evaluated.  Xt  is  designed  to  synthesise  the  spectrum  of  any  low- 
frequency  even  function,  but  it  could  easily  be  modified  to  synthesise 
the  spectra  of  odd  functions.  The  principles  are  also  applicable  to  a 
design  which  could  synthesise  the  spectrum  of  a  general  function.  In 
the  latter  case,  either  a  two -channel  memory  or  a  pair  of  memories 
is  required,  and  the  synthesiser  becomes  considerably  more  complex. 

The  Fourier  synthesiser  operates  upon  signals  with  frequencies 
which  range  from  dc  to  an  upper  limit  that  is  imposed  by  the  allowable 
sampling  rate.  The  present  upper  limit  is  0.  5  cps.  An  improvement  by 
a  factor  of  five  is  expected  as  soon  as  several  components  in  the  system 
are  replaced. 

The  system  has  been  designed  for  use  in  Fourier  spectroscopy  to 
synthesise  the  spectral  distribution  from  the  corresponding  interferogram. 
It  has  been  tested  in  thiB  application,  and  the  performance  is  satisfactory. 
It  was  used  with  a  lamellar  grating  interferometer,  but  it  could  also  be 
U6ed  with  a  well-compensated  Michelson  interferometer. 


J. 
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Additional  work  is  required  to  reduce  the  errors  to  a  minimum. 

This  may  be  accomplished  by  adding  an  apodizaticm  system  and  by  changing 
several  of  the  components  in  the  synthesizer.  A  more  accurate  multiplier, 
a  better  analog  to  digital  converter,  and  a  larger  memory  with  more 
addresses  will  improve  the  accuracy.  A  system  to  interpolate  the  data 
which  is  in  the  memory  will  also  improve  the  performance. 

A  secondary  design  criterion  has  been  met.  This  was  to  design  a 
system  that  would  be  portable  and  one  that  could  also  be  used  outside 
a  laboratory.  The  present  synthesizer  can  meet  this  requirement  by 
being  packaged  into  a  more  compact  form.  The  choice  of  the  type  of 
summation  unit  and  the  decision  to  use  a  crystal -controlled  frequency 
synthesizer  and  special  gating  circuitry  instead  of  a  phase  locked 
oscillator  chain  permitted  the  construction  of  a  fairly  compact  system. 

The  exclusive  use  of  salid  state  circuits  also  made  this  small  size 


possible 
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APPENDIX  1 

Design  Equations  and  Data  for  the  Nonlinear  Shaping  Networks 


The  90-degree  detector  used  several  non-linear  networks  to 
shape  a  siaewavc  into  a  series  of  narrow  pulses.  Biased  diode  networks 
are  used,  and  the  components  are  chosen  to  give  a  network  with  an 
approximate  quartic  or  quin  tic  transfer  characteristic.  Thus, 


i  =  k  E4 
out  in 


or 


i  =  kE. 
out  in 


where  k  is  chosen  to  give  the  components  reasonable  values.  Six  line 
segments  are  used  to  approximate  the  characteristic,  and  the  breakpoints 
are  chosen  at  e  =  1,  3,  5,  7#  and  9  volts.  Figure  39  shows  part  of  the 
characteristic  and  the  first  few  line  segments  of  the  approximation. 

Let  the  voltages  at  the  breakpoints  be  E^,  E^,  E^»  E^,  E^,  and  E, , 
respectively.  Conductances  will  be  used  in  the  equations  as  it  simplifier 
the  calculations.  The  values  will  then  be  converted  to  resistances, 
which  are  shown  in  Figure  23  .  At  the  first  breakpoint,  E^,  the 
current  for  the  quintic  network  is 


When  e  =  E^,  it  is 


5 

i  =  G  E,  =  k  E, 
o  1  1 


IOC 


i  *  goe2  +  g1(e2-e1)  =  k  E2 

Likewise,  at  e  * 

i  *  goe3  +  G1{E3-E1)  +  G2(E3-E2)  =  k  e3 
Similarly,  at  e  - 

i  *  goe4  +  g1{e4-e1)  +  g2{e4-e2>  +  g3{e4-e3>  =  k  E4 

Also,  at  e  =  E_ 

5 

i  r  GoE5  +  Gj(E5«Ej)  +  G2(E5-E2)  +  G3(E5-E3) 

+  g4(e5-e4)  =  k  e4 

Finally,  at  e  * 

i  =  gqe6  +  g1(e6-e1)  +  g2(e6-e2)  +  g3(e,-e3) 

+  G4<E6'E4>  +  G5(E6-E5>  *  kE6 

The  equations  are  solved  in  turn,  and  the  conductances  (resistances)  a 

determined.  The  value  of  k  is  chosen  such  that  the  largest  resistance 
7 

is  less  than  10  ohms. 
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A  basic  diagram  of  a  non-linear  network  is  shown  in  Figure  40 
while  Figure  41  shows  a  practical  network  that  has  been  used*  The 
breakpoint  voltages  are  determined  with  the  aid  of  a  resistance  voltage 
divider*  One  percent  precision  resistors  are  employed  and  zener  diodes 
are  used  to  control  the  voltages  at  the  ends  of  the  divider*  The  resistors 
in  the  divider  have  low  values,  compared  to  the  ones  in  the  shaping 
network.  This  is  necessary  if  the  breakpoint  voltages  are  to  remain  nearly 
constant  as  the  input  voltage  varies. 


APPENDIX  2 


110 


Analysis  and  Performance  of  Two  Preliminary  Methods 
for  Obtaining  Fourier  Transforms 

Electro-Optical  Spectrum  Analyzer 

An  electro-optical  computer,  somewhat  similar  to  Montgomery* s  f !  0 j , 
was  constructed*  In  this  case,  an  oscilloscope  wao  used  to  furnish  both 
the  light  source  and  the  harmonic  functions.  Montgomery  used  a  projector 
lamp  and  sine-cosine  masks  of  different  frequencies. 

A  periodic  function  can  be  expanded  in  a  Fourier  series  of  the  forri 


f(t)  = 


00 

E 

n=l 


(a  cosu  t  +  b  sinw  t) 
'  n  n  n  n 


where 


a 

n 


T 


T/2 

f  m  costo  t  dt 
-T/2 


and 


T/2 

b  e  ~  r  f(t)  sinu  t  dt 

n  7  '  w  n 

-T/2 

For  a  non-periodic  function,  the  Fourier  transform  is 


Ot 

F(u)  =  f  F{t)e'jutdt 


Ill 


The  functions  to  be  examined  are  either  even  or  odd.  Then  only  one  type 
of  terms,  cosine  or  sine,  is  needed  in  the  series  expansion.  For  non¬ 
periodic  functions,  only  the  cosine  or  the  sine  transform  is  required. 

The  multiplication  and  integration  may  be  performed  optically.  In 
the  following  diagram,  light  entering  from  the  left  is  modified  by¬ 
passing  through  a  mask  whose  transmission  function  is  T.  The  light 
between  the  masks  can  be  represented  as  I^T^  while  to  the  right  of  the 
second  mask  it  is  I  T.T,.  Thus,  the  light  output  is  the  product  of  the 
two  transmission  functions. 
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The  transmission  functions  are  masks  of  variable  density  or 
variable  width.  Two  variable  density,  or  one  variable  density,  and 
one  variable  width  function  may  be  multiplied. 

If  the  transmission  function  of  one  of  the  masks  is  f(x)  and  the  other 
mask  carries  the  sine  or  cosine  function,  the  integral  is  obtained  by 
focusing  the  emerging  light  onto  a  photocell.  The  output  of  the  photoce’l 
is  a  constant  times  the  integral 


f(x)  cosuxdx  or  /  £{x)  sinwxdx. 
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Several  Fourier  coefficients  may  be  obtained  by  the  use  of  several  slides, 
each  with,  a  cosine  or  sine  function  of  a  different  "spatial"  frequency. 

The  experimental  system  is  shown  in  Figure  42a  .  The  100  kc 
oscillator,  operating  in  conjunction  with  a  slow  sweep  speed,  produces 
a  band  of  light  on  the  screen*  The  vertical  gain  is  adjusted  so  only  the 
most  uniform  central  region  of  the  pattern  is  used.  The  variable  frequency 
oseillator  (VFO)  is  connected  to  a  motor  which  drives  the  tuning  mech¬ 
anism.  The  frequency  is  varied  from  100  cps  to  20  kc  in  a  period  of 
several  minutes.  This  oscillator  is  used  to  intensity  modulate  the  set  pe 
in  a  sinusoidal  manner.  A  trigger  circuit  produces  synchronizing  pul  ies 
to  start  successive  sweeps  at  the  same  point  with  respect  to  the  intensify 
modulating  signal.  As  the  frequency  increases,  fringes  of  increasing 
spatial  frequency  move  across  the  screen. 

The  mask  carries  the  £(x)  function.  The  light  emerging  from  th<. 
mask  is  focused  onto  a  photocell.  The  cell  output  is  amplified  end 
recorded  on  a  strip  chart  recorder. 

Several  test  functions  were  evaluated  <o  determine  the  system 
performance.  Variable  width  and  variable  density  functions  were  used., 
Variable  width  masks  were  made  by  modulating  a  carrier  with  the  der.irc 
function  and  photographing  the  amplitude  modulated  wave  from  a  scope, 
using  Polaroid  transparency  film.  Variable  density  masks  were  wade 
by  projecting  a  Ronchi  Ruling  onto  photographic  film  with  an  enlarger 
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and  developing  the  film*  A  sample  of  the  analysis  of  the  Ronchi  Ruling 
data  is  shown  in  Figure  43a  ,  and  a  sample  of  the  variable  width  data  is 
shown  in  Figure  43b  . 

This  system  has  many  disadvantages  which  render  it  unfeasible  for 
an  accurate  high  speed  Fourier  analysis  technique: 

1.  The  number  of  fringes  that  can  be  placed  on  the  cathode  ray 
tube  is  limited  by  resolution  and  defocuaing  effects.  With  100  fringes 
on  the  screen,  defocusing  at  the  edges  is  very  noticeable  while  the 
central  region  is  still  in  good  focus. 

2.  The  intensity  is  not  uniform  and  varies  across  the  face  of  the 
scope.  Thij  is  called  the  "shading"  characteristic. 

3.  A  very  high  noise  level  existB  due  to  the  continuously  varying 
pattern.  As  each  new  fringe  comes  on  the  screen,  the  total  light  inteu. 
varies.  This  is  especially  evident  when  only  a  few  fringes  are  present 
and  causes  very  large  variations  in  the  detected  signal. 

4.  Higher  harmonics  in  the  spectrum  are  lost  as  a  result  of  the 
combined  effects  of  noise  and  defocusing. 

5.  Because  of  the  high  noise  level,  amplitude  comparison  between 
different  frequency  components  is  difficult, 

6.  The  intensity  modulation  is  not  uniform  as  a  function  of  frequc  c 
and  good  cosinusoidal  intensity  distributions  are  difficult  to  ofct?in. 
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7.  Two  closely -spaced  frequencies  in  a  width  modulated  pattern 
are  not  accurately  resolved, 

8#  Initial  phasing  of  the  scope  pattern  and  the  mask  is  done  by  eye 
and  is  not  sufficiently  accurate.  Errors  cause  the  sine  function,  repre¬ 
senting  a  harmonic,  to  be  asymmetrical. 

9,  The  method  is  not  suited  for  the  real-time  analysis  of  an 
electrical  signal. 

A  Photographic  Spectrum  Synthesizer 

Experiments  were  performed  to  determine  the  feasibility  of  using 
photographic  techniques  to  compute  Fourier  transforms,  A  diagram  cr 
the  system  is  ahown  in  Figure  42b  ,  The  VFO  is  manually  controlled, 
and  a  counter  is  used  to  determine  the  discrete  frequencies  required  in 
the  synthesis.  A  set  of  polarizers  controls  the  intensity  of  the  light 
which  reaches  the  film.  The  camera  photographs  each  succeeding 
pattern  as  a  multiple  exposure  on  the  same  frame. 

The  function  under  investigation  is  even  and  is  represented  in 
Figure  44  ,  It  is  properly  sampled,  and  the  sampled  function  is  alee 
shown*  The  sampled  function  may  be  represented  as 
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It  is  possible  to  obtain  the  Fourier  transform  of  F(t)  by  -uperposi  :on. 
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where  Fn(u>)  is  the  Fourier  transform  of  F^{t).  Conversely 


Fo(u)  +  F1(u)  +  F2<w)  +  ..Fn(u)  =  f  Fo{t)e*‘Jwtdt  + 


+  J  Fn(t)e  a  dt 


Thus 


F(u)  a  £  Fq(u) 
n=0 

The  function  has  been  reduced,  by  sampling,  to  a  set  of  delta 
functions,  and  the  total  transform  is  a  .■  animation  of  the  transforms  c  f 
the  delta  functions. 

The  Fourier  transform  of  a  6  function  at  t  =  0  is  of  constant 
amplitude  and  extends  to  oo  in  both  directions  in  the  frequency  domain. 
The  Fourier  transform  of  a  pair  of  delta  functions,  symmetrically  placer’, 
about  t  =  0,  is  a  cosine  wave  whose  frequency  depends  on  the  location  of 
the  6  functions  from  t  =  0. 
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The  function  to  be  analyzed  is  manually  sampled,  and  the  sample 

amplitudes  are  recorded*  These  are  then  converted  into  the  angles  at 

which  the  crossed  polarizers  must  be  set  when  the  photographic  data  are 

taken.  These  control  the  intensity  of  the  light  reaching  the  film,  and 

each  sample  corresponds  to  a  different  intensity  and  spatial  frequency. 

The  film  is  developed,  dried,  and  processed  with  a  micro-densitometer 

which  plots  the  density  of  the  negative  versus  position.  The  more  dense 

portions  are  points  where  the  fringes  have  added  in  phase  and  represent 

harmonics  of  the  original  function. 

Several  simple  functions  were  analyzed*  Data  are  presented  in 

Figures  45  and  46  which  show  the  results  of  analyzing  a  i - - — ! 

*  IT  t  j 

pulse,  a  finite  cosine  wave  and  a  simple  interferogram.  The  system  has 
disadvantages  which  render  it  unacceptable  for  a  real-time  synthesizer: 

1*  The  samples  are  manually  obtained  and  converted  to  polarizer 
angles. 

2.  The  system  is  feasible  for  only  even  or  odd  functions, 

3.  Film  grain  adds  greatly  to  the  noise  level. 

4.  It  is  difficult  to  maintain  the  shape  of  the  cosinusoidal  fringes 
as  the  frequency  is  altered, 

5.  The  pattern  does  not  have  uniform  focus. 

6.  There  is  a  non-uniform  shading  characteristic  that  cannot  be 
exactly  compensated. 
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7 .  The  film  is  non-linear  and  tends  to  Saturate  as  a  large  number 
of  samples  are  taken* 

8.  The  time  involved  between  the  acquisition  of  the  original  data 
and  the  spectrum  is  excessive. 

9.  The  accurate  control  of  the  polarizer  setting  is  difficult. 

Because  of  the  large  number  of  disadvantages,  it  was  decided  to 

abandon  both  schemes  as  a  method  of  obtaining  Fourier  transforms  in 
real  time.  An  investigation  was  made  to  determine  the  feasibility  of 
using  storage  tubes  as  a  summation  device,  but  they  have  all  of  the 
disadvantages  of  a  CRT  and  are,  therefore,  not  acceptable. 

The  system  that  is  used  is  an  all -electronic  analog  of  the  photo¬ 
graphic  system.  It  has  the  advantages  of  being  automatic,  of  utilizing  a 
linear  summation  scheme  that  will  not  saturate,  and  it  is  able  to  genera 
pure  cosine  functions  of  the  desired  frequency.  The  electronic  scheme 
is  also  a  real-time  technique. 
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